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SUMMARY 

Two-wavelength  interferometry,  employing  either  a  ruby  (6943-2)  or  a 
neodymium  (10,600a)  laser  with  a  second  harmonic  generator,  was  used  to  study 
argon  plasmas  in  a  shock  tube  for  shock  Mach  numbers  from  17  to  24 ,  initial 
pressures  from  3.0  to  0.3  mm  Hg  and  for  degrees  of  ionization  up  to  50$.  The 
measured  equilibrium  electron  densities  for  both  wavelength  pairs  were  about 
1.9$  greater  than  the  calculated  values  accounting  for  the  excited  states,  but 
1.6$  less  than  the  values  based  on  the  neglect  of  these  states.  Assuming  a 
pla/.ma  of  free  electrons  and  ground-state  atoms  and  ions,  the  interferograms 
yielded  ratios  of  the  polarizability  for  the  argon  ion  to  that,  for  the  atom  of 
■0.650  +  .033  at  3471a  and  0.647  +  .04l  at  5300a.  These  results  are  in  general 
agreement  with  the  value  of  0.716  for  the  ratio  of  the  static  polarizabilities 
obtained  by  the  method  of  Slater  screening  constants.  The  uniformity  of  the 
measured  electron  densities  and  polarizabilities  with  change  in  wavelength 
suggests  that  the  effect  of  the  excited  state  refractivities  are  not  important 
for  temperatures  up  to  13,500°K. 
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1.1  Introduction 

A  knowledge  of  the  electric  dipole  polarizability  for  atomic  of  molecular 
species  is  an  essential  prerequisite  to  obtaining  a  better  understanding  of  how 
individual  atoms  and  molecules  interact  both  with  ea,ch  other,  with  free  electrons 
and  other  forms,  of  matter  and  finally  with,  static  and  dynamic  electric  and  mag¬ 
netic  fields.  Typical  properties  which  are  intimately  related  to  the  polariz¬ 
ability  are  the  long  range  intermolecular  dispersion  force  constant,  the  Verdet 
constant  of  Faraday  optical  rotation,  the  diamagnetic  susceptibility,  the  Rayleigh 
scattering  cross  section,  the  total  effective  oscillator  strength  and  the  die¬ 
lectric  constant. 

However,  it  is  the  relationship  between  the  polarizability  a  for 
species  s  and  the  phase  refractive  index,  p  representing  the  bulk  property  of 
the  transmitting,  medium  which  is  of  primary  interest  In  the  present  circumstance; 
for  dilute  gases  this  relationship  is  given  by  (Ref.  3) 

(Ug  -  1)  =  27T  jsag  (1.1.1) 

where  N  is  the  particle  number  density  of  species  s. 

S  '  *  '  ;  > 

Experimentally  measured  refractive  indices  provide  a.  relatively  accurate 
yardstick  against  which  the  merits  of  the  presently  available  theoretical  cal¬ 
culations  can  be  judged.  However,  a  more  practical  need  arises  during,  the 
optical  and  partictiiafiy  the  interferometric  study  of  transparent  systems  which 
are  undergoing  either  physical  or  chemical  changes,,  clearly  the  refractive 
indices  of  all  species  present  during  the  reaction  must  be  known- prior,  to  inter¬ 
pretation  of  the  experimental  data  .obtained  through  interferometry. 

In  particular,,  this  laboratory  is  presently  undertaking  an  inter¬ 
ferometric  study  of  nonequilibrium  shccx  generated  argon  plasmas  which  requir  es 
a  knowledge  of  the  refractivities  of  all  the  plasma  species  present,  viz.*,  the 
free  electrons,  the  neutral  argon  atoms  and  the  singly  charged  ions.  As  the. 
refractivity  for  this  latter  species  is  .as  yet  unknown,  it  is  to  be  the  purpose 
of  the  present  study  to  experimentally  determine  a  value  for  this  quantity., 

1.2  Review 

Previous  measurements  of  the  polarizabilities  of  chemically  active 
gaseous  species,  normally  existing  only  at  high  temperatures,  have  been  made 
almost  exclusively  using  the  shock  tube  in  conjunction  with  a  (Mach-Zehndor) 
interferometer  (Refs.  1  to  9).j  the  authors  and  respective  species  are  indicated 
in  Table.  1. 

Although  numerous  optical  and  particularly  interferometric  studies 
have  been  made  on  argon  plasmas,  the,  contribution  of  the  bound  .electrons  to  the 
overall  plasma  refractive  index  is  usually  neglected  or  at,  best  approximated  . 
by  the  neutral  atom  ground  state  value  in  view  of  the  relatively  large  re¬ 
fractivity  and  strong  dispersion  of  the  free  electrons . 

There  have  been  two  studies  which  have  specifically  set  out  to 
measure  the  refractivity  of  an  argon  plasma.  Firstly,  Alpher  and  White  (Refs. 1,2), 
who  pioneered  the  use  of  the  shock  tube-interferometer  technique'  in  measuring 
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the  polarizability  of  high  temperature  species,  employed  interferometry  simul¬ 
taneously  at  two  wavelengths.  By  assuming  the  dispersion  for  the -atom  and  ion 
to  be  negligibly  small,  they  showed  that  the  measured  and  calculated  electron 
densities  were  in  good  agreement,  but  with  the  average  experimental  value- being 
about  2 %  greater  than  the  value  calculated  from  the  shock  Hugoniot  equation. 
From  this  agreement,  it  was  concluded  that  the  theoretical  free  electron  re¬ 
fractive  index  given  by  (see  Section  2.i), 

(p  -  1)  =  -  4.484  x  10" l2*  A2  N  ,(1*2.1) 

e  e 

where  A  and  N  are  the  wavelength  and  electron  particle  density  respectively, 
is  a  good  description  pf  reality  in  the  optical  region  of  the  spectrum*.  The 
reason  for  this  small  but  systematic  discrepancy  between  the  measured  and' 
calculated  electron  densities  is  not  clear.  It  appears  that  their  theoretical 
calculation  employed  Saha  and  energy  equations  using  only  the  ground  state 
terms  in  the  summation  for  the  electronic  partition  functions  of  the  neutral 
atom  and  first  ion.  Consequently  this  calculated  value  is  an  upper  bound  (see 
Section  4.3)  suggesting  that  the  discrepancy  steins  from  some  other  unspecified 
source.  In  addition,  Alpher  and  White  indicated  that  the  measured  piasina 
refractive  index  was  consistently  larger  than  the  calculated  value  but  gave 
no  numerical  estimate  for  this  difference.  This  calculation  assumed  a.  plasma 
of  free  electrons  and  ground- state  atoms  and  ions  where  the  ion  refractivity 
was  taken  to  be  0.67  of  the  corresponding  neutral  atom  value  as  obtained 
using  Slater  screening  constant  theory**  (Ref.  10,).  However,  it  should  be 
noted  that  this  theory  is  at  best  only  approximate  as  indicated  by  the  fact 
that  the  absolute  value  of  the  neutral  argon  ground  state  static  polarizability 
is  about  2|  times  the  well  documented  experimental  value.  Alpher  and  'White 
suggested  several  possible  explanations  for  this  observed  difference  between 
the  measured  and  calculated  plasma  refractive  indices  and  in'  .particular  that 
the  excited  states  might  be  contributing  to  the  experimental  measurements. 

It  is  therefore  concluded  from  the  work  of  Alpher  and  White  that  the 
figure  of  O.67  for  the  ratio  of  the  polarizability  of  the  first  ion  to  that 
for  the  neutral  atom  is  of '  the  correct  order  of  magnitude,  but  that  due  tp  un¬ 
certainties  ip  the  measured  plasma  refractive  index  and  electron  density,  it  Is 
not  possible  to  place  a  figure  on  the  accuracy  of  this  result. 

The  second  study  was  made  by  Hug,  Evans,  Tankin  and.  .Camibel  (Refs. 

13,  14)  who  employed  optical  interferometry  to  study  axisymmetric  argon  plasmas 
produced  by  either  a  plasma  jet  or  a  free  burning  arc.  Experimental  values 
for  the  plasma  refractive  index  were,  found  to  be  generally  about  12%  below  the 
calculated  value  based  on  a  plasma  of  free  electrons  and  ground  state  atoms 
and  ions,  it  was.  suggested  that  this  might  be  due  to  the  effect  of  the  highly 
excited  bound  electrons  lying  close  to  the  continuum  behaving,  as  if  free  there¬ 
by  producing  a  reduction  in  the  value  of  the  observed  plasma  refractive  index. 
However,  it  should  be  mentioned  that  Alpher  and  White  made  their  measurements 

*  Alpher  and  White  (Ref.  2)  gave’  a  value  for  the  .constant  in  Eq... 1.2.1  equal 
to  4.46  whereas  the  present  author  calculates  a  value >of  4*484  using  the 
self  consistent  set  of  physical  constants  given  in  Ref .. -78.  ■ 

y  y 

Both  the  present  author  and  Ascpli-Bartoli  et  al  (Ref.  11*)'  calculate 
a  value  of  0.72  for  this  ratio  using  the  same  theory  as  given.iin  Ref. 

10.  In  addition  Alpher  and  White  quote  a  value  of  O.69  in  Ref.  12 
using  the  same  theoretical  model. 
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at  4 122$  and  5463$  whereas  Hug  et  al  used  a  He-Ne  laser  at  6328$  (Eel's'.1. 13,  l4) 
and  also  at  11,523$  (Ref.  15)  so  that  the  latter  measurements  were  more  sensitive 
to  the  free  electron  concentration.  It  would  therefore  seem  that  the  experiments 
of  Hug  et  al  were  not  sufficiently  sensitive  to  resolve  the  contribution  of  the 
argon  ion  to  the  total  plasma  refractive  index. 

It  is  therefore  the  purpose  of  the  present  work  to  measure  the  induced 
electric  dipole  polarizability  for  singly  ionized  argon  by  employing  short  wave¬ 
length  interferometry  to  study  a  shock  tube  generated  argon  plasma. 

The  shock  tube  is  capable  of  providing  a  known  concentration  of  a.  given 
species  at  a  well  defined  thermodynamic  state  over  a  wide  range  of  temperature 
and  pressure.  In  addition,  the  gas  sample  is  one-dimens iDnal  such  that  reduction 
of  the  interferometric  data  is  most  straightforward. 

The  refractive  indices  of  the  atom  and  ion  plasma  species  are  virtually 
independent  of  wavelength  whereas  the  value  for  the  free  electrons  given  by  Eq. 
1.2.1,  is  strongly  dispersive  dictating  that  the  experimental  measurements  be 
made  at  as  short  a  wavelength  as  is  consistent,  with  the  experimental  system; 
this  ideal  is  readily  achieved  by  using  the  second,  harmonic  of  the  ruby  laser 
whose  wavelength  of  3^71$  lies  just  below  the  cut-off  frequency  for  conventional 
glass  optics. 

In  addition,  by  paying  careful  attention  to  accuracy  in  the  calculation 
and  solution  of  the  shock  Hugoniot  equations  and  in  the  experimental  measurements, 
the  ion  refractivity  was  measured  with  a  respectable  degree  of  accuracy  in  spite 
of  the  dominant  contributions  of  the  free  electron  and  neutral  atom  terms  to 
the  overall  plasma  refractivity. 

Finally,  both  theoretical  and  experimental  approaches  are  made  to  de¬ 
termine  whether  the  neutral  atom  excited  states  contribute  significantly  to  the 
measured  plasma  refractive  index. 

2.  THEORETICAL  CONSIDERATIONS 

2.1  Plasma  Refractive  Index 

For  pressures  of  the  order  of  one  atmosphere,  the  plasma  refractive 
index  can  be  given  by  (Ref.  12) 

(Up-i)p,T  =  ^s'^PjT  (2^1.1) 

s 

where  the  summation  extends  over  all  ’s'  species  which  constitute  the  plasma  at 
pressure  p  and  temperature  T.  The  assumption- of,  additivity  implies  that  the  - 
individual  plasma  components  do  not  interact  with  each  other,-  viz,  that,  the 
long  range  attractive  intermolecular  fprces  do  npt  contribute  to  the  refractivity^ 
Both  experimental  and  theoretical  values  of  the  polarizability  for  argon  at 
high  pressure  show  that  deviations  from- .the  zero  pressure  value  arenno  , greater 
than  l/2$  at  100  atmospheres  (Ref.  16).  It  is  therefore  assumed  that  the  re¬ 
fractivity  for  a  given  species  exhibits  a  linear  dependence  ofi  density  (Refi.  17). 

For  a  plasma  of  neutral  atoms,  first  ions  and  free  electrons,  Eq. 

2.1.1  can  be  expanded  to  give 


(2.1.2) 


^e_1^p,T 


^X-l)p,T 


where  the  free  electron,  neutral  atom  and  first  ion  terms  are  denoted  by  sub¬ 
scripts  'e',  'A'  and  'I*  respectively.  The  neutral  atom  and  first  ion  components 
implicitly  contain  the  contributions  fi:om  the  excited  electronic  states  in 
addition  to  that  for  the  ground  state. 


Ditchburn  (Ref.  l8)  has  given -the  classical  derivation  for  the  re¬ 
fractive  index  p  of  a  gaseous  species,  .based  on  the  dipole  moment  induced  in 
the  bound  optical  electron  oscillator  by  an.  indident  electromagnetic  wave; 
this  is  given  by  (Ref.  l8)  1 


(•*!  - 1)  - 


f 

r 
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where  f  is  the  number  of  electron  oscillators,  in -each  atom  with  resonant  fre¬ 
quency  e  is  the  electron  charge,  m  is  the  electron  mass,.  Ng  is  the  number 
density  of  species  s,  v  is  ttje.  frequency  of  the  incident  electrpmagnetic  .radia¬ 
tion,  and  y  ^  is  the  absorption  damping  factor  for  electron  oscillator  type  r . 


In  the  optical  region  of  the  spectrum  two  sinplifications  can  be  made 
to  Eq.  2.1.3.  Firstly,  y is  very  small  compared  to  y so  that  provided  v, 

the  frequency  of  the  incident  radiation,  is  remote  from  the  resonance  frequency 
v  ,  then  the  absorption  damping  term  can  be  omitted.  Secondly,  in  the  optical 
region  p  for  gaseous  species  generally  has  a  value  close  to  unity  so  that 
Eq.  2.1.§  can  be  reliably  approximated  as 
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(2.1.4) 


For  a  free  electron  gas,  the  resonance  frequencies  are  all  zero  so 
that  Eq.  2.1.4  becomes 

-e8N  p 

(p  -l)  =  - -  =  -  4,4844  x  10-14  N  A2  (2.1.5) 

2my  e 

where  N  is  the  electron  density  and  A  (=c/v)  is  the  wavelength  of  the  incident 
radiation.  Alpher  and  White  (Refs.  1,  2)  who  measured  the  free  electron  re- 
fractivity  in  the  optical  region  found  their  results  to  be  in  close  agreement 
with  Eq.  2.1.5  (see  Section  1.2). 


The  quantum  mechanical  analogue  of  Eq,  (2.1.4)  for  the  species  re¬ 
fractive  index  is  essentially  the  same  except  that  the  bound  -electron  is  now 
a  quantized  oscillator  making  transitions  between  any  two  states;®  and  n, 
rather  than  a  fixed  oscillator  r,  Due  to  the.  existence  of  an  infinite  array  of 
possible  quantum  states,  the  quantized  analogue  of  Eq.  .2.1. ‘4  'must 'be.  obtained 
by  summing  qver  all  m  states.  This  formulation'  was  first  given  by  Ladenberg 
(Refs.  19,20)  who  showed  that  for  dilute  gases  in  the  absence-  of  external 
fields,  the  total  refractivity  for  a  system  of  electrons  bound  in  atomic  and 
ionic  configurations  for  incident  radiation  of  frequency-  v,  can  be  given  by 
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(2.1.6) 


where  m  and  n  refer  to  any  two  discrete  energy  levels,  m  being,  the  lower  of  the 

two  states  whose  statistical  weights  are  g„  and  respectively.  N  and  Ef- 

m  n  •  m  n 

are  the  number  densities  of  the  atoms  in  these  two  states  and  v  ,  is.  the  frequency 

,  '  mn 

of  the  radiation  emitted  in  the  transition  from  n  to  m  for  which  the  absorption 
oscillator  strength  is  f^.  'similarly,  subscripts  L  and  c  refer  to  energy  levels 
existing  at  the  ionization  limit  and  in  the  continuum  respectively. 

The  second  term  in  the  expression  for  the  discrete  transitions  repre¬ 
sents  the  effect  of  negative  dispersion  or  induced  emission  which  only  becomes 
important  at  high  tenperatures  where  the  populations  of  the  excited  states  be¬ 
come  significant.  It  should  be  noted  that  the  experimental  demonstration  of 
the  existence  of  this  phenomenon  by  Ladenberg  (Ref.  19)  was'  an  important  confir¬ 
mation  of  the  correctness  of  the  theory  of  quantum  mechanics . 

The  integral  within  the  curly  brackets  in  Eq.  2.1.6  represents  the 
contribution  to  the  species  refractive  index  due  to  transitions'  from  discrete 
states  M  to  continuum  states  C,  lying  at  energies,  greater  than  hvQ^  above  the 

ground  state  where  hv^  is  the  energy  required  to  ionize  a  neutral  ground  state 
atom.  The  partial  derivative  within  the  integral  is  related  to-  the  phot.oipni- 

zation  cross-section  cr  (v  )  associated  with  absorption,  from  discrete  state  M 

m  me  .  -  • 

to  continuum  state  G  by  the  expression  (Ref.  23), 


=  ov(,  )  as, 

■dv  m'  me'  -  2 

me  ire 

It  should  be  noted  that  these  two  phenomena  are  both  predicted  by  the  quantum 
theory  of  matter  and  have  no  analogue  in.  the  classical  theory. 


At  temperatures  such  that  the  populations  of  the  excited  states 
are  small,  Eq.  2.1.6  can  be  simplified  so  that  the  species  refractive  index 


p  can  be  given  by  the  ground  state  term  for  which  m  =  0,  viz 
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where  N  is  the  ground  spate  number  density  and  is  hence  equivalent  to  the  -total 
s 

species  density  in  the  low  temperature  limit. 

Using  the  few  available  resonance  oscillator  strength  values  for 
neutral  argon  (Ref.  21),  the  refractivity  for  this  species  was  calculated 
using  the  first  term  of  Eq.  2.1.7  and  was  found  to  constitute  only  2 k°jo  of  the 
experimental  value  (Ref.  22)  at  3^71^,  the  wavelength , of  the  second  harmonic 


I 


of  the  ruby  laser.  Similarly,  the  corresponding  value  calculated  for  singly 
ionized  argon  is  only  9$  of  the  value  calculated  from  Slater  screening  constant 
theory  (see  Section  1.2)  where  (|ij-l)/(n^-.l).  S  O.67.  Clearly  many  more  resonance 
oscillator  strengths  for  discrete  transitions  must  be  known  before- this  approach 
can  be  used  to  obtain  the  refractive  indices  of  species  having  relatively  complex 
spectra  such  as  argon.  In  addition,  resonance  transitions  to  the  continuum  states 
can  be  expected  to  make  a  significant  contribution  as  is  the  case  for  helium 
(Ref.  23).  As  yet  no  wave  mechanical  calculation  appears  to  have  been  made  for 
t.ie  dispersion  of  either  neutral  or  singly  ionized  argon. 

With  regard  to  the  experimental  determination  of  the  dispersion  for 
neutral  argon,  the  situation  is  more  satisfactory.  Two  recent  measurements, 
the  first  by  Peck  and  Fisher  (Ref.  22)  for  visible  and  near  infrared  wavelengths 
and  the  second  by  Chaschina,  Gladushchak  and  Schr eider  (Ref.  24)  for  vacuum 
pltra- violet-  wavelengths ,  both  show  very  close  agreement  when  the  data  of  Peck 
and  Fisher  is  extrapolated  into  the  vacuum  ultra-violet  region  (Ref.  24).  Peck 
and  Fisher  gave  their  data  in  the  form  .of  the  empirical  dispersion  formula 

(nA-i)6  =  6.786711  x  icf5  +  giP^2  .  (2.1.8) 

A  0  (1.1*4  x  ICX^A2-  1). 

where  the  subscript  'O'  denotes  conditions,  of  standard  density  and  temperature 
(STP) ^  viz.  0°C  and  760  mm  Hg  and  where  A  is  the  wavelength  of  the  incident 
wavelength  in  Angstrom  units  ( R ) .  As  these  measurements  were  made  for  ambient 
conditions,,  the  populations  of  the  excited  states  are  therefore  insignificant, 
so  that  Eq.  (2.1.8)  ip  effect  represents  the.  refractivity  for  .neutral  ground 
state  argon  and  is,  in  principle,  equivalent  to  the  expression  given  by  Eq. 
(2.117).  The  experimental,  data  represented  by  Eq.  (2.1.8)  will  therefore  be 
used  to  describe  the  dispersion  of  neutral  (ground  state)  argon  as  it  exists 
in  the  preshock  state  for  laser  wavelengths  from  the  near  infra-red  through 
to  the  near  ultra-violet  regions  of  the  spectrum  (see  Sections  3  and  4). 

However,  before  Eq.  2.1.8  can  be  used  to  describe  the  neutral  atom 
refractivity  component  of  the  plasma  refractive,  indei  in  Eq.  2.1.2,  it  has  to 
be  shown  that  the  refractivity  for  neutral  argon  remains,  for  all  practical 
purposes,  constant  for  temperatures  up  to  13,500°K  typical  of  the  experiments 
to  be  described  in  Section  4.  This  implies  that  the  effects  of  both  negative 
dispersion  and  the  excited  state  polarizabilities  must  make  an  insignificant 
contribution  to  the  overall  species  refractivity  in  the  wavelength  region  of 
interest. 


In  principle  the  species  refractivity  given  by  Eq.  2.1.6  is  a  function 
of  temperature  due  both  to  the  strong  temperature  dependence  of  the  population 
number  densities  N  and  to  the  fact  that  the  excited  state  polarizabilities 
and  statistical  weights  diverge  rapidly  with  increasing  principal  quantum, number. 


Firstly,  consideration  is  given  to  tb?  contribution  of  negative  dis¬ 
persion  to  the  overall  species  refractivity.  If  the  atomic  and  ionic  electronic 
energy  levels  are  assumed  to  exist  in  a  Boltzmann  distribution,  then  the  popu¬ 
lation  number  densities  are  given  by 


Ng  gffi  exp  (-E/KE) 

CO 

Y  giexp("VkT') 


(2.1.9) 


V 


6 


where  the  summation  in  the  denominator  is  the  electronic  partition  function  for 
species  's',  Ng  is  the  total  number  density  for  species  's'  and  E  in  the  energy 
of  the  mth  electronic  state  above  the  ground  state  (m  =  0) . 

For  the  ground  state,  the  negative  dispersion  factor  in  Eq.  2.1.6  is 
(l-Nn/NQ.  go/gn)  which  assumes  the  value  (i  rexp(-En/kT))  for  a  Boltzmann  distri¬ 
bution.  A  'worst  case'  estimate  of  this  term  can  now  be  obtained  by  .giving  Effi 
a  constant  value  corresponding  to  the  lowest  excited  state  which  for  neutral  , 
argon  is  that  for  the  3Pp  metastable  level  (3p^4s  configuration)  lying  931^3  cm”1 
above  the  ground  state  (Ref.  25).  This  term  is  now  a  constant  with  respect  to 
m  and  can  therefore  be  taken  outside  the  summation  in  Eq.  2.1.$;  at  13,500°K  this 
factor  has  a  value  of  (1  -  2  x  10"5)  which  is  essentially  equal  to  unity. 

As  the  contribution  of  the  excited  state  polarizabilities  to  the  overall 
species  refractivity  is  itself  considered  to  be  of  second  order  at  least  for  the 
range  of  temperatures  for  the  present  experiments,  it  is  felt  that  the  overall 
effect  of  negative  dispersion  on  these  values  will  be  small.  However.,  as  the 
energies  involved  in  discrete  nonrresonance  transitions  fof  neutral  argon  can  be 
considerably  smaller  than  those  for  the  resonance  transitions,  it  was  fblt  pru¬ 
dent  to  include  this  correction  term  into  the  calculation  made  to  estimate  the 
contribution  of  the  excited  state  polarizabilities  to  the  neutral  argon  refrac¬ 
tivity.  This  calculation  will  be  dealt  with  shortly. 

An  indication  of  the  magnitude  of  polarizability  for  the  excited  states 
can  be  gauged  from  the  value  for  the  first  excited  state  of  neutral  argon,  i.e., 
the  3P2  level  of  the  3p^s  configuration.  Bederson  et  al  (Refs.  26,  27)  have 

used  an  electric  and  magnetic  field  balance  technique  to  measure  the  polarizability 
for  this  metastable  state  and  found  it  to  have;  a  value  about  30  times  greater  than 
that  for  the  ground  state.  Unfortunately  this  technique  is  only  suitable  for 
measuring  the  polarizabilities  of  metastable  excited  stateg  due  to  their  rela¬ 
tively  long  lifetimes.  It  should  be  noted'  thatj  at  13,500  K,  the'  population  of 
the  3Pg  level  is  about  1/5000  of  the  ground  state  value  so  that  considered  alone, 
this  contribution  is  negligible.  However,  for  the  hydr ogeni 6 - typ e  gas,  the  polari¬ 
zability  for  the  q^h  principal  quantum  level  varies  as  q°  .(Ref.,  5)  •  Although 
this  functional  relationship  is-  not  exactly  true  for  the  case  'of  argon,,  it  does 
demonstrate  the  very  strong  dependence  of  the.  polarizability  on  "the  principal 
quantum  number.-  This  point  receives  added  emphasis  when  it  is  -realized  that  the 
q  =  21  level  with  a  statistical  weight  g  =  3200  may  be  occupied-  for  an.  argon 
plasma  at  13,500°  K  and  1  atmosphere  pressure  prior  to  termination  of  the  series 
by  the  plasma  microfield  (Ref.  28).  It  should,  however,  be  noted  that  although 
the  plasma  microfield  restricts  the  occupation  of  the  discrete  levels  to  -a.  finite, 
number,  it  is  pow  possible  to  have  transitions  to  continuum  states  at  energies 
which  formerly  existed  below  the  undepressed  ionization  limit..  As  a  near  cori-r 
tinuous  array  of  discrete  terminal  levels  has  now  been  replaced  by  a  continuous 
array,  it  is  suggested  that  the  effect  of  the  .plasma  microfield  on  the  species 
refractive  index  is  negligible. 

An  attempt  has  therefore  been  made  to  estimate  the-  summed  effect,  of 
the  excited  state  polarizabilities  for  neutral  argon  by  applying-  Eq.  2,1.6  Vith. 
m  =  1  through  to  infinity  together  with  the  available  tabulations  of  o.scillator 
strength  values  for  discrete  non-resonance  transitions.  A  total  of  657  terms 
was  obtained  from  the  recent  National  Bureau  of  Standards  compilation  (Ref..  21) 
of  most  likely  values,  together  with  additional  theoretical  values  from  the 
calculations  of  Murphy  (Ref.  29)  and  of  Johnston. (Ref .30)  and  additional  experi¬ 
mental  values  from  the  experiments 
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of  Wiese  et  al  (Ref.  79)* 

Comprehensive  calculations  of  the  polarizability  for  any  one  excited 
electronic  state  using  this  method  requires  a  detailed  knowledge  of  a  very  large' 
number  of  oscillator  strength  values ,  including  the  contribution  from  discrete- 
continuous  '  ransitions .  The  present  calculation  can  therefore  at  best  only  hope 
to  estimate  the  dispersive  effect  due  to  the  stronger  discretertransitions  in  the 
neutral  argon  spectrum. 

The  contribution^ of  discrete  non-resonance  transitions  to  the  total 
excited  polarizability  contribution  is  obtained  from  Eq.  2.1.6  in  t.ie  form 
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•For  a  Boltzmann  distribution  of  level  populations,  this  equation  can  be  rewritten 
as: 


(h  *  -1)  = 


00  CO 


v  A  “1=1  n=m+l  '  ■  mn 


(2.1.11) 


where  A  corresponds  to  transition  n.  to  m,  N.  is  the  neutral  argon  atom  particle 

ion  ^  A  . 

density, 

f  g  Q  2  2 

A 
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in  the  Einstein  coefficient  of  spontaneous  emission  and-  is  the  partition  function 
for  neutral  argon. 

The  calculation  was  performed  for  a  number  of  ‘.air  wavelengths  in  the  . 
near  ultra  violet, visible  and  near  infrared'  regions  and  in.  particular  for  those 
wavelengths  corresponding  to  the  ruby  laser,  A  =  6943*5  8  (Ref.  .80  )  and  its  second 
harmonic,  to  the  neodymium  laser,  A  =  10,600A‘  and  its  second  .harmonic  and  to  the 
He-Ne  laser  wavelengths  at  63288  and  11523 A. 

Whereas  the  spectral  bandwidths  for  the  ruby  arid.,  heliumnneon  lasers  are 
generally  less  than  0.1a,  that  for  the  neodymium  laser  corisists  of  a  system  of 
lines  distributed  over  a  50A  bandwidth.  Unfortunately  the  spectral  structure  of 
the  neodymium  laser  output  is  a  function  of  a  number  of  laser  parameters,  e.g. , 
pump  energy,  Q-switch  mode,  laser  glass  type  arid' quality,  cavity  end  reflectors, 
so  that  it  is  difficult  to  control  or  predict  the  spectral  nature  of  the  neodymium 
laser  radiation.  However,  there  appears  to  be  as  yet  no- oscillator  strength 
data  available  corresponding  to  any  dominant  neutral  argon  transitions  which  fall 
in  the  spectral  regions  of  the  neodymium  laser  or  its  second  harmonic.  Hence  the 
use  of  an  average  value  for  the  wavelengths  of  the  neodymium  laser  and  its  second 
harmonic  is  felt  to  be  permissible. 

For  wavelerigths  shorter  than  10,0008''  '-.Wright  et.  al  (Ref..  31)  have 
tabulated  both  the  allowed  and  the  forbidden  transitions  for  neutral  argon.. 

From  these  tables  it  is  apparent  that  the  only  allowed  transition  lying  close 
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to  any  of  the  visible  6r  ultra-violet  laser  lines  is  that  in  the  5d-4s  array 
with  an  air  wavelength  of  3^ 72. 54.8.  Unfortunately  the  oscillator  strength  for 
this  transition  has  been  neither  calculated  nor  measured.  However,.  Desai  and 
Corcoran,  who  measured  oscillator  strengths  for  neutral  argon  in  this  spectral 
region  (Ref.  32),  have  noted  that  the  emission  for  this  line  lies  at  or  below 
the  noise  level  of  their  measurements  such  that  the  Einstein  spontaneous  emission 
coefficient  for  this  transition  has  a  value  £  1.3  x  10“5  sec"l  (Ref.  33).  This 
value  was  therefore  included  in  the  calculation  for  the  total  excited  state 
refractivity  although  a  preliminary  calculation  indicated  that  its  contribution 
at  3471. 75S  was  not  large. 

Although  there  are  a  large  number  of  neutral  argon  transitions  in  the 
region  of  the  He-Ne  laser  line  at  II523S,  no  similar  tabulation  exists  for  the 
transitions  in  the  vicinity  of  this  laser  line.  The  results  for  this  wavelength 
should  therefore  only  be  considered  as  tentative  although  none  of  the  Known 
transitions  lie  close  to  this  wavelength. 

The  calculation  was  performed  for  temperatures  of  10,0Q0°K,  I3>500°K, 
15,000  K  and  20,000°K  for  a  plasma  pressure  of  one  atmosphere  with  the  corres¬ 
ponding  neutral  atom  partition  function  taken  from  the  tabulation  of  Drellishak 
et  al  (Ref.  28). 

The  total  excited  state,  refractive,  index  contributions  for  the  above 
mentioned  laser  wavelengths  and  plasma  conditions  are  shown  in  Table  2  where,  the 
major  contribution  occurs  at  6943,5a. 

For  piasma  conditions  of  13,50O°K  and  1  atmosphere  pressure  . at  6943.5$, 
the  refractivity  for  the  excited  electronic  states  of  neutral  argon  is  about 
6.6$  of  the  total  species  value  or  about  6.1$  of  the  corresponding  ground  s;tate 
value .  . 

Figure  1  shows  the  total  neutpal  argon  refractivity  as  a  function  of 
wavelength  for  the  stated  experimental  conditions  which  in  addition  to  the -ground 
state  contributions,  includes  those  either  from  the  total'  of  6.57  non-resonance, 
discrete  transitions  or  from  the  39  transitions  of  the  4p-4s  array,  ,  The  ground 
state  contribution  to  the  overall  neutral  atom  refractivity  as  measured  by  Peck 
and  Fisher  (Ref.  22)  is  also  shown  for  comparison. 

„  It  immediately  becomes  obvious  that  for  wavelengths  less  than  about  . 
9000A,  the  transitions  in  the  4p-4s  array  make  the  dominant  contribution  to 
the  excited  state  refractivity  and  are  responsible  for  the  broad  resonance  extrema 
which  are  centered  at  about  8200A  and  9100A.  This  is  not  surprising  as  the 
transitions  for  this  array  tend  tp  fall  into  two  wavelengthsgroups  centered  on 
the  transitions  having  vacuum  wavelengths  of  8117.5$  and  9125-8  which  have  by 
far  the  largest  (g^)  values  in  their  respective  groups. 

For  longer  wavelengths,  the  influence  of  a  great  many  relatively  strong 
near  infrared  transitions  becomes  apparent  by  vittue  of  the  irregular  nature 
of  the  data  shown  in  Fig.  1.  However,  the  fact  that  the  shorter  wavelength, 
data  appears  to  be  continuous  is  purely  a  result  of  choosing  the  diagnostic 
wavelengths  fo  be  sufficiently  remote  from  any  individual  transitions.,  i.e.,, 
(?-?')>  20a  say,,  so  that  local  resonances  do  not  ocpur.  (Equation  2,1.11  is 
by  its  very  nature  a  non-continuous  function  of  the  wavelength  of  the  incident 
radiation) „ 

An  interesting  observation  regarding  the  two  extrema  shown  in 
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Fig.  1,  is  that  they  are  not  due  to  a  single  transition  but  to  about  20  out  of 
the  total  of  30  transitions  in  the  4™-4s  array.  In  addition,  these  transitions 
all  originate  from  the  resonance  levels  which  for  an  equilibrium  plasma  are  the 
most  highly  populated  of  the  excited  states  of  neutral  argon. 

It  should  be  possible  to  confirm  the  existence  of  these  resonances  by 
performing  single  frequency  interferometry  on  an  equilibrium  plasma  of  known 
thermodynamic  and  ground  state  refraetivity  properties  at  a  number  of  different 
wavelengths  in  the  neighbourhood  of  these  extrema.  An  ideal  source  of  diagnostic 
radiation  for  such  measurements  could  be  provided  by  a  pulsed  dye  laser  which 
can  be  tuned  to  operate  in  this  wavelength  region. 

If  the  existence  of  these  exbrema.is  confirmed,  then  it  should  be 
possible  to  monitor  directly  the  population  density  of  the  rate-limiting  resonance 
states  for  neutral  argon.  This  could  be  achieved  in  a  manner  similar  to  that 
proposed  by  Measures  (Ref.  34) ,  who  suggests  that  by  performing  dual  frequency 
interferometry  at  wavelengths  close  to  and  relatively  remote  from  the  line  cehter 
for  a  single  transition,  one  can  separate  out  the  fringe  shift  contributions  due 
to  changes  in  the  other  unknown  plasma  species-  from  that  due  to  the  single  tran¬ 
sition  under  -observation.  In  the  present  situation,  due  to  the  broad  spectral 
nature  of  these  extrema,  the  fringe  shifts  at  both  wavelengths  would  contain 
contributions  due  to  the  transitions  in  the  4p-4s  array  as  well  as  those  due  to 
the  free  electrons  and  the  ground  st.ates  of  neutral  and  singly  ionized  argon. 
However,  this  is  no  obstacle  to  interpretation  of  the  data  as  only  those  com¬ 
ponents  due  to  the  excited  state  transitions  wili  exhibit  a  marked  wavelength 
dependence  over  a  relatively  small  wavelength  interval.  For  example,  at 
8000$  with  AA  =  4o$,  the  change  in  the  free,  electron  refraetivity  is  about  T/o, 
which  to  a  first  approximation  can  be  ignored. 

Although  these  calculations  appear  to  indicate  that  the  excited  states 
make  a  small  but  real  contribution  to  the  neutral,  argon  refraetivity  for  tempera¬ 
tures  greater  than  10,000  K  and  particularly  for  the  ruby-  laser  fundamental 
wavelength  of  6943. 5a,  it  was  decided,  to  omit  this  excited  state  contribution 
from  the  neutral  argon  refraetivity.  'This  was  due  to  the  uncertainty  introduced 
into  the  calculation  by  the  omission  of  a  large  number  of  unknown  oscillator 
strength  values  for  transitions  between  bound  states  in  addition  to  those  for 
discrete-continuum  transitions.  Rather,  it  was  decided  to  use  the  experimental 
data  of  Peck  and  Fisher  (Ref.  22)  given  by  Eq.  2.1.8,  to  describe  the  total 
neutral  argon  refraetivity  for  the  shock  generated  argon  plasma  at  pressure  p 
and  temperature  T. 

As  this  excited  state  refraetivity  contribution  should  selectively 
effect  the  experimental  determination  of  the  free,  electron  density  as  obtained, 
from  the  two-wavelength  interferometric  technique  (see  Section  2.3)*  which  is 
primarily  dependent  on  the  fringe  shift  at  6§43.5$  and  not  on  the  fringe  shift 
at  3471 .75^,  then  a  comparison  between  the  measured,  electron  density  made  using 
the  ruby  and  neodymium  lasers  and  the  corresponding  calculated  values  should 
reflect  the  presence  of  a  significant  excited  state  refraetivity  contribution. 
This  comparison  will  be  discussed  in  Sections  4.3  and  4.4’  in  the  light,  of  the 
present  experimental  measurements. 

In  view  of  the  relatively  small  neutral  argon  excited  state  refrac- 
tivity  contributions  to  the  plasma  refractive  index  for  the  chosen  diagnostic- 
wavelengths,  a  similar  calculation  to  estimate  the  effect  of  the  excited  states 
of  singly  ionized  argon  does  not  seem  to  be  necessary.  The  population  of  the 
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first  excited  state  for  this  species  at  13,500.  K  is  approximately  3  x  10"  of  the 
ground  state  value  which  is  considerably  smaller  than  for  the  resonance  states 
of  neutral  argon.  In  addition,  a  major  obstacle  to  such  a  calculation  exists  in 
that  no  oscillator  strength  data  is  available  for  the-  non-resonanb  transitions 
from  this  level,  although  numerous  values  are  available  for  transitions  between 
the  more  energetic  but  sparsely  populated  levels  (Ref.  21). 


The  expression  for  the  plasma  refractive  index  to  be  used  in  the  analysis 
of  the  present  experimental  measurements  will  therefore'  be  given  by  Eg.  2.1.? 
with  the  free  electron  term  provided  by  Eq.  2.1.5,  the  neutral  atom  term  By .Eg 12 $1. 8 
and  with  the  ion  refractivity  term  to  be  determined  as  described  in  Sections  2.2 
and  b.4.  ' 


2.2  Relationship  Between  the  Fringe  Shift  and  the  Change  in  the. Plasma 
Phase  Refractive  Index 

The  fringe  shift  S^  as  obtained  using  a  Mach-Zehnder  interferometer  is 
related  to  the  change  in  the  phase  refractive  indices  between  any  two  states  1  and 
2,  by  the  relationship. 


(2.2.1) 


where  superscript  A  refers  to  the  wavelengths  of  the  diagnostic  radiation  and  L 
is  the  geometrical  path  length  over  which  the  change  occurs.  Expanding  the  sum¬ 
mation  terms  in  Eq.  2.2.1  into  the  form  given  by  Eq.  2.1.2  arid  denoting  the  pre¬ 
shock  and  post  shock  (plasma)  states  by  subscripts  1  and  2  respectively,  the 
fringe  shift  equation-  now  becomes  ' 


(  x)'[  t(|1A-1,2t(Mr1>2-<:(,X€-1)2  {(PA-I)l}l 


(2.2. 2) 


where  it  is  assumed  that  the  preshock  gas  at  state  1  consists  -entirely  of  neutral, 
argon  atoms . 

As  the  component  refractivites  are-  linearly  dependent  on  their  res¬ 
pective  densities?,  then- Eq.  2,2.2  can.be  rewritten  in  terms  of  the  refractivities 
at  standard,  temperature  and  .pressure  -,(STP)  to  give- 


( x )-( 

N  n  *  >• 


(a,2*3) 


where  subscript  ’o’  denotes  conditions  at  STP,  x^is  the  degree  of  single  ioni¬ 
zation,  p  is  the  mass  density  and  (p^-l)o  and  (p^l)Q  are  given  by  Eqs.  2.1.8 

and  2.1.5  respectively. 

The  expression  to  be  used  to  determine  (ji^-l)  ,  the  refractivity  for 
singly  ionized  argon  at  STP,  is. obtained  directly  from  Eq.  2.2.3as 
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where  x^, 


the  degree  of  ionization  is  eliminated  through  use  of.  the  equality 


iq  _o 

and  Ng  =  2.687  x  10  7  cm  is  Loschmidt's  number.  The  quantities  to  be  experi¬ 
mentally  determined  are  therefore  Sl'p,  the  fringe  shift  between  states  1  and  2, 
p  and  p0  the  preshock  and  equilibrium  plasma  mass  densities  and  N  „  the  plasma 

free  electron  density. at  stab:  2. 

In  the  optical  region  of  the  spectrum, the  refractivities  for  the  heavy 
particles  (neutral, atoms,  ions)  are  effectively  non-dispersive  whereas  the  free 
electron  refractivity  varies  as  A2-;-  .  It  is  therefore  'dear  from  Eq..  2.2.3  that, 
the  fringe  shift  contribution  for  the  heavy  particles  varies  as  A“i  whereas:  that 
for  the  free  electrons  varies  as  A.  Hence  in  order  to  increase  the  fringe  shift 
contribution  due  to  the  argon  ion  and  consequently  the  sensitivity  of  the  experi¬ 
mental  measurements  while  simultaneously  depressing  the  contribution  due  to  the 
free  electrons,  it  is  essential  to  conduct  the  experiment  at  as  short  a  wave¬ 
length  as  possible.  As  will  be  described  in  Sections  3  and  b,  this  is  accompli¬ 
shed  by  employing  the  second  harmonic  of  the  ruby  laser  at  3^71* 75a  which  lies 
just  below  the  cut-off  frequency  for  conventional  glass  optics.  ' 


2.3  Determination  of  the  Plasma  Electron-' Density  by  the  Method  of. Two 
wavelength  Interferometry 


In  order  to  obtain  (pT-l) 

A  .  I  o 

the  fringe  shift  ,  the  values  of 


from  Eq,  we  must  know  in  addition  to 

P2  and  Xg.  These  can  be  obtained  from,  the- 


shock  Hpgoniot  equations  for  given  values  of  the  shock  Mach  number  Mi  ,  p  and 
T  on  the  assumption  that  the*  gas  is  in  Saha  equilibrium.  However,  in  tne 
presence  of  radiation  cooling,  it  cannot  be  categorically  stated  that  .the  shock 
generated  plasma  attains  Saha  equilibrium  (Refs.  58,  59>  60).  Clearly  Saha 
equilibrium  will  exist  if  the -electron  density  calculated  via  the  shock 
Hugoniot  equations  is  in  agreement  with  the  corresponding  experimentally 
measured  value.  It  is  the  determination- of  this  Tatter  value  with  wtf.;ch  we 
are  now  concerned. 


As  discussed  in  Section.  1.2  and  2.1 j  the  two  wavelength  interferometric 
technique  of  Alpher  and  White  (Refs.  '2,  12)  is  capable  of  providing  an  accurate 
measure  of  the  electron -density  for  both  equilibrium  and  non  equilibrium  plasmas. 

Two  simultaneous  interferograms.  of  the  same  shocked  gas.  sample  are 

obtained  for  wavelengths  A  and  A  as  far  apart  as  possible,  with  A  >  A  at  as 

x  y  y  x 

long  a  wavelength  as  the  plasma  or  optical  system  allows.  The  1;  4">er  require¬ 
ment  reflects  the  strong-  dispersive  nature,  of  the  free  electron  refractivity 
where  (p^-l)  a  A2  in  contrast  to-  the  weak  dispersion  of  the  heavy  particle  r.e- 

f ractivities .  By  simultaneously  solving  Eq.  2.2.3  for  wavelengths  A  -and  A 
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where  it  is  explicitly  assumed  that  the  excited  states  do  not  contribute  to  either 
the  neutral  atom  or  the  ion  refractivities-..  The  term  in  curly  brackets  in-  the 
denominator  of  Eq.  2.3.1  represents  the  difference  between  the  ratio  of  the  ion 
refractivity  to  the  neutral  atom  refractivity  at  A^  to  that  at  A^„  The  resonance 

absorption  transitions  for  neutral  argon  occur  at  1048.2$  and  1066. 7^  (Ref.  2l) 
with  the  result  that  the  dispersion  for  neutral  argon  is  very:  weak,  in  the  near 
ultra-violet,  visible  and  near  infrared  regions  of  the  spectrum.  The  correspon¬ 
ding  transitions  for  singly  ionized  argon  occur  at  918. 88  and  932. oX  (Ref;  21) 
so  *hat  the  dispersion  for  this  species  can  be  expected  to  be  even  weaker  than 
that  for  neutral  argon  in  the  same  spectral  region.  Hence  this,  difference  term 
referred  to  above  is  likely  to  be  very  small.  However,  by  taking  the  pessimisti¬ 
cally  large  value  pf 

7=  ±  °*1  -  ; 
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which  for  the  wavelengths  of  the  ruby  laser  (69433$)  arid,.-  (3471-.7$)respond  to  about  a 
i6$  change  in  the  ion  refractivity,  the  denominator  in  Eq.  2.3.1  takes  a  value 
of  (l  +  0.0C63).  Similarly  for  the  wavelengths  of  the  neodymium  laser  (10600a) 
and  its  second'  harmonic  ( 5300A) ,  this  term  has  a  value  of  (l  +  Q.0027),.  It  is 
therefore  clear  that  this  term  has  a  negligible  effect  oh  the  electron  density 
given  by  Eq,  2.3.1  under  normal  conditions  and  can  hence  be  omitted. 


The  expression  to  be  used  to  obtain  the  experimental  electron  density 
from  the  fringe  shift  data  is  therefore  given  by 
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Two  points  regarding  Eq.  2.3.2  call  for  comment.  Firstly  this  expression 
is  completely  independent  of  the  (as  yet  unknown);  ion  refractivity  and  only  weakly 
dependent  on  thennehtral  atom  refractivity  through  the  ratio 

>  ■ 

which  has  a  value  close  to  unity  in  the  optical  region.  By  assuming  this  ratio 
to  be  unity,  an  error  of  0.8%  is  generated  in  We  for  the  wavelengths  of  the  ruby, 
laser  and  its  second  harmonic.  For  the  neodymium  laser  and  its  second  harmonic, 
the  corresponding  error  of  0.2$  is  less  significant. 
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Secondly,  the  accuracy  with  which  Ne  cari  be  measured  is  dependent  on  the 

magnitude  of  (A  -A  )/(A  +  A  ) .  For  a  laser  (A  )  and  its  second  harmonic  (A  ) 
y  x"  y  x  y  x 

this  terms  takes  a  value  of  l/3 ♦  As  (A  -A  )  approaches  zero,  then  both  the  numera- 

y  x 

tor  and  denominator  approach  zero  so  that  the  value  Ne  becomes  highly  susceptible 
to  errors  in  and  .  It  should  be  mentioned  that  Alpher  and  White  (Ref.  2, 
12)  made  the  additional  assumption  that  (p^y-l)  =  (pV^-l)  =  const,  with  the 

result  that  0  ° 

1,8 "  7VT  (i?)[4-4(^}][(^)  -1]  (2-3-3) 


"y  .  \  "  /  L  ^  \  y  /  J'  L  \  V  /  J 
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However,  in  their  case,  the  wavelength  separation  was  somewhat  smaller 
with  (A  -  A^)/(A  +  Ax)  =  0.l4,  so  that  this  approximation  is  more  acceptable 

although  a  price  is  paid  in  accuracy  as  has  already  been  noted. 


3.  EXPERIMENTAL  FACILITIES 


3.1  Shock  Tube 

A  detailed  account  of  the  construction  and  operation  of  the  UTIAS  4"'  x 
7"  combustion  driven  shock  tube  has  been  given  by  Boyer  (Ref.  35) •  However, 
several  modifications  and  improvements  were  made  to  this  facility  with  a  view  to 
expediting  the  present  experiments.  These-  are  described  in  the  following  sections. 

3-1.1  Shock  Tube  Driver  Section 

-  -  -  -  i  -  4 

A  number  of  refinements  have  been  made  in  the  operation  of  the  combus¬ 
tion  driver  facility.  It  was  noticed  that  the  stpichiometric  mixture  of  hydrogen 
and  oxygen  diluted  by  75$  of  helium  tended  to.  exhibit  rough  combustion.  Thisi 
became  manifest  as  a  series  of  high  frequency  (~  lKhz),  moderate  amplitude  or.cii- 
lations  superimposed  on  the  main  pressure  rise,.  The  reason  for  the  pressure  u.races 
given  by  Boyer  (Ref.  35)  not  exhibiting,  these  oscillations, .may  be  related  to  the  . 
change  in  length  in  the  driver  tube  from  the.  original  value  of  13.75  ft.  to  the 
present  length  of  50  in.  It  should  be  mentioned  that  the  period  of  these  pressure 
pulses  corresponds  closely  to  that,  of  a  single,  sound  wave- reflecting  between  the 
end  walls  of  the  50  in.  long  combustion  chamber  .cavity  at  the  sound  speed  of  the 
75$  helium  combustion  mixture  (Ref.  36) • 

As  it  is  conceivable  that  these  pressure  pulses  might  amplify  into,  a 
detonation  wave  or  alternatively  damage  and  possibly  break  the  ignition  wire,, 
a  change  to  80$  helium  dilution  was  made  in  view  of  the  experiences  of  Nagamatsu 
and  Martin  (Ref.  37).  Although  this  measure  had  the  desired  effect  of  achieving 
smooth  combustion,  it  was  also  observed  to  result  in  a,  small  reduction  in  shock 
Mach  number  for  otherwise  identical  initial  conditions  together  with  a  much 
increased  time  to  peak  combustion  pressure,  viz,  100  msec  as  against  10  msec. 

There  are  several  possible  explanations  for  this  reduction  in  per¬ 
formance.  Firstly,  the  more  rapid  rise  time  at  75$  helium  dilution  may  tend  to 
accelerate  the  diaphragm  opening  process  thereby  promoting  more  efficient  shock 
formation.  Secondly,  a  longer  residence  time  fob  combustion,  in  principle* 
leads  to  higher  thermal  losses  with  a  consequent  drop  in  pressure  and  sound,  speed 
of  the  driver  gas.  Finally  and  significantly,  calculations  (Ref.  36)  indicate 
that  the  peak  combustion  pressure  at  80$  helium  dilution  is  about  10$  less  than 


for  75 %  dilution  for  the  same  predomination  pressure-. 

A  compromise  dilution  of  77g$  was  therefore  tried  and  found  to  give 
smooth  combustion  for  which  the  risetime  to  peak  combustion,  pressure  was  com¬ 
parable  to  that  for  7%  helium  dilution.  Although  Benoit  (Ref.  36)  indicates 
that  this  degree  of  dilution  achieves  a  peak  pressure  approximately  3lo  less  than 
that  for  the  7%  dilution,  a  similar  loss-  of  shock  tube  performance  was  not 
observed.  %  , 


In  addition,  with  a  view  to  ensuring  complete  combustion,  the  7lM° 
dilution  was  partitioned  into  components  of  %  hydrogen  and  72g$  helium.  This 
ensures  better  hydrogen- oxygen  contact  and  more  rapid  combustion  while  simul¬ 
taneously  raising  the  post-combustion  sound  speed  without  creating  any  undue 
hazard  due  to  the  presence  of  an  unburnt  hydrogen  fraction  in  the  dump  tank. 

t  ' 

For  this  77§$  level  of  h-lium  dilution,  the  experimentally  determined 
shock  tube  characteristic  for  combustion  driving  into  argon  using-  a  single  'dia¬ 
phragm  was  found  to  be  given  by  in  (p^/pj)  -  0.422  Mg  +  4v02)  for-  10:  <  Mg  <  24., 

where  p^  is  the  peak  combustion  pressure,  p^  is-  the  preshock  (initial)  test  gas 
pressure  and  is  the  incident  shock  Ma,ch  number.  The  ratio  of  p^  to  pw.y,  the 

pre- combustion  pressure  for  the  above  conditions,  was  found-  to  hayenan 
experimental  value  of  6.85  +  0.25  whereas  the  calculated-  value  was  about  8.05 
(Ref.  36).  As  p^  is  initially  determined  by  the  shock  requirements,  the  con¬ 
servative  yalue  of  (Pi|/Pmj_x)  =  6.60  was  therefore  employed  in  conducting  the 
present  experiments. 

The  other  change  made  to  the  combustion  system  involved  increasing  the 
thickness  of  the  tungsten  heating  wire  used-  to  ignite  the  combustipn  .mixture  from 
a  diameter  of  0.010"  to  0.015"  with  a  view  to  reducing  the  possibility -of  a  .break 
before,  during  or  after  a  run  and  hence  to  ensure  use  of  the  wire  for  seyer.al 
runs.  Initially,  the  50  inch  length  of  0.015  in.dia.  wire  has  a  resistance,  of 
0.6l  ohms  which  tends  po  double  rafter  each  run  due  to  oxidation.  The  wire  was 
usually  changed  after  the  resistance  exceeded  2  ohms  or  after  -every  three  runs 
whichever  comes  first,.  For  high  press  we  runs  (p^>  5000  psi)  however,  the  wire 
tended  to  break  during  combustion  but  without  any  otherwise  deleterious  effects 
on  the  combustion  process  or  the  facility. 

In  order  to  heat  the  tungsten  wire  to  the  same  temperature  in  order 
to  maintain  smooth  combustion,  the  voltage  on  the  4§pF  capacitor  was  increased 
by  50$  to  12.75  Kv  providing  an  energy  density  of  about  90  Joules  per  foot  of 
0.015"  diameter  wire,  sufficient  to  make  the  wire  glow  a  bright  red-orange  colour 
when  tested  in  air  at  one  atmosphere  pressure. 

A  common  cause  of  detonation  is  localised  ignition.  In  the  present  , 
system  these  conditions  were  observed  to  exist  in  the  form  of  arcing  at  the 
wire  terminals  due  to  inadequate  clamping  pressure.  This  became  particularly 
important  in  view  of  the  above  mentioned  voltage  increase.  Careful  redesign 
of  the  end  clamps  was  found  to  be  sufficient  to  eliminate  this  localised  spark¬ 
ing. 

3.1.2  Criteria  for  Design  of  the  Shock  Tube  Diaphragms 

The  bursting  pressure  characteristics  for  circular  scribed  diaphragms 
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given  by  Boyer  (Ref.  35)  refer  to  discs  of  annealed  stainless  steel  type  302  with 
an  unsupported  diameter  of  6|  inches.  However,  the  diaphragms  now  used  have  been 
standardized  to  stainless  steel  type  £04  with  a  2B  finish  for  which  the  diaphragm 
clamp  has  a  clear  diameter  of  7-1/8  inches.  The  production!  of  new  calibration 
curves  was  therefore  found  to  be  necessary. 


A  series  of  tests  using  the  hydraulic  diaphragm  bursting  rig  described 
in  Ref.  35  was  performed  on  annealed  diaphragms  of  type  304  stainless  steelhhaving 
thicknesses  of  0.062,  0.109  and  0.172  inches.  The  experimentally  measured  burst¬ 
ing  pressures  for  the  scribed  diaphragms  p  was  adequately  described  by  the 
dimensionless  relationship  ^ 

(|M0n 

where  p^  is  the  bursting  pressure  of  the  unscribed  diaphragm,  h  is  the  residual 
diaphragm  thickness  after  scribing  and  t  is  the  total  diaphragm  thickness.  The 
measurements  indicated  a  value  for  n  =  2.20  +  0.1  for  (h/t)  >  0.10  where  the 
error  limits  for  n  are  due  to  a  number  of  uncertainties,  e.g.,  variations  in 
material  composition,  degree  of  anneal,  diaphragm  wedge  angle  and,  total  thickness 
h  and  also  to  errors  in  measurements  of  p^,  p^  and  h. 

■2- 

The  bursting  pressure  of  the  unscribed  diaphragms,  p^.,  was  found  go  be 
given  by  the  usual  relationship  (Ref.  81) 


(3.1.2 .2) 


where  cr  is  the  ultimate  tensile  strength  for  the  annealed  material(.85000  psi  for 
SS304),  d  is  the  unsupported  diameter  of  the  diaphragms  and -K  is  a  factor  which, 
corrects  for  the  non- ideality  of  the  diaphragm  performance  due  to  the  competing 
effects  of  work  hardening  and  diaphragm  thinning  which  occur  during  deformation 
just  prior  to  bursting.  Ideally  K  should  be  unity  and  in  reality  !ha3,’-'a  value  close 
to  this  value.  Unfortunately  a  universal’  experimental  value  of  K  was-  not  ob¬ 
tained,  possibly  due  to  variations  in  material  composition,  degree  of  anneal, 
wedge  angle  and  surface  finish.  Specific  values  of  K  for  each  diaphragm  thick¬ 
ness  are  as  follows: 


t  =  0.172  in,  K  =  1.01 

t  =  0.109  in,  K  =  1.10  j-  for  SS304  and  d  =  7-1/8.  inn 

t  =  0.062  in,  K  =  1.06  *  ' 


For  combining  Eq.  3*1*2.1  and  3*1. 2. 2  the  bursting  pressure  for  a  given  scribed 
diaphragm  is  given  by  •  .  .  . 


(3.1.2j3) 


which  when  combined  with  the  experimental  values  of  n  and  K,  gives  a  reliable 
prediction  of  diaphragm  performance  .for  (h/t;)  >  0.10. 
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an  unsupported  diameter  of  6|  inches.  However,  the  diaphragms  now  used  have  been 
standardized  to  stainless  steel  type  -304  with  a  2B  finish  for  which  the  diaphragm 
clamp  has  a  clear  diameter  of  7-1/8  inches.  The  production  of  new  calibration 
curves  was  therefore  found  to  be  necessary. 


A  series  of  tests  using  the  hydraulic  diaphragm  bursting  rig  described 
in  Ref.  35  was  performed  on  annealed  diaphragms  of  type  304  stainless  steelhhaving 
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where  p^  is  the  bursting  pressure  of  the  unscribed  diaphragm,  h  is  the  residual 
diaphragm  thickness  after  scribing  and  t  is  the  total  diaphragm  thickness.  The 
measurements  indicated  a  value  for  n  =  2.20  +  0.1  for  (h/t.)  >  0.10-  where  the. 
error  limits  for  n  are  due  to  a  number  of  uncertainties,  e.g.,,  variations  in 
material  composition,  degree  of  anneal,  diaphragm  wedge  angle  and,  total  thickness 
h  and  also  to  errors  in  measurements  of  p^,  p^  and  h. 

The  bursting  pressure  of  the  unscribed  diaphragms,  p^,  was  found  to  be 
given  by  the  usual  relationship  (Ref.  81) 
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where  e  is  the  ultimate  tensile  strength  Cor  the  annealed  material(85000  psi.  for. 
SS304),  d  is  the  unsupported  diameter  of  the  diaphragms  and  K  is  a  factor  which 
corrects  for  the  non- ideality  of  the  diaphragm  performance  due  to  the  competing 
effects  of  work  hardening  and  diaphragm  thinning  which  occur  during  deformation 
just  prior  to  bursting.  Ideally  K  should  be  unity  and  in  reality  ‘has  a  value  close 
to  this  value.  Unfortunately  a  universal  experimental  value  of  K  was  not  ob¬ 
tained,  possibly  due  to  variations  in  material  composition,  degree  of  anneal, 
wedge  angle  and  surface  finish.  Specific  values  of  K  for  each  diaphragm  thickr 
ness  are  as  follows: 
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For  combining  Eq.  3.1»2.1  and  3*1*2.2  the  bursting  pressure  for  a  given  scribed 
diaphragm  is  given  by  '  .  . 


t  =  0.172  in,  K  =  1.01 

t  =  0.109  in,  K  =  1.10  J  for  SS304  and  d  = 

t  =  0.062  in,  K  =  li06  •  ' 


(3. 1.2:3) 


which  when  combined  with  the  experimental  values  of  n  and  K,  gives  a  reliable 
prediction  of  diaphragm  performance. for  (-h/t;)  >  0.10. 


3.1.3  Shock  Tube  Driven  Section 


Modifications  to  the  driven  section  were  primarily  centered  around  im¬ 
proving  the  shock  tube  vacuum  with  a  view  to  obtaining  a  high  purity  test  gas 
sanple  in  the  test  section  prior  to  shock  arrival. 

Through  the  introduction  of  a  6  in.  dia.  vacuum  diffusion  pump  between 
the  shock  tube  and  Roots  blower  vacuum  pump  (see  Ref.  35) >  a  vacuum  and  leak- 
outgassing  rate  of  about  1  x  10”5  mmHg.  and  1  x  10“5  mm  Hg/min.  respectively  were 
obtained  after  a  24  hour  pumpdown  period.  The  relatively  long  pumpdown  period 
is  due  to  the  large  volume  of  the  shock  tube  dump  tank  (35*4  cu.ft.)  in  relation 
to  the  shock  tube  channel  (9.6  cu.ft.)  together  with  the  low  conductance  of  the 
shock  tube  which  connects  the  dump  tank  to  the  vacuum  pumps . 

If  desired  the  vacuum  and  combined,  outgassing-leak  rate  could  be  im¬ 
proved  by  about  an  order  of  magnitude  by  extending  the  pumpdown  period  -from 
2  to  3  days.  However,  as  events  transpired,  a  high  degree  of  tes.t-gas  purity 
was  not  required  particularly  as  the  argon  test  gas  was  intentionally  seeded 
with  approximately  0.4$  of  molecular  hydrogen  by  pressure  for1  reasons  to-be 
discussed  in  Section  4.2. 

The  argon  and  hydrogen  test,  gases  used  inotheae  experiments  -were  both, 
of  the  'Linde  High  Purity  Grade'  type.  The  argon,  has  a  manufacturers  minimum 
purity  of  99*996$  by  volume  with  a  typical  trace  gas  analysis  of  less  than  15 
ppm  of  nitrogen,  7  ppm  of  oxygen  and  5  ppm  of  carbonaceous  matter  together  with 
a  dew  point  of  better  than  -71°F.  Similarly  the  hydrogen  has  a  quoted  purity  of 
99*99$  by  volume  with  less  than  100  ppm  of  nitrogen  and  10  p*  l  of  oxygen  to¬ 
gether  with  a  dew  point  better  than^-71°F. 

The  hydrogen  was  first  admitted  into  a  calibrated  reservoir  to  a  pre¬ 
determined  pressure  and  then  released  into  the  shock  tube.  The  argon  test  gas 
was  then  added  directly  to  the  shock  tube  until  the  desired  total  initial 
pressure  was  attained. 

3*1*4  Measurement  of  the  Incident  Shock  Wave  Properties 

The  measurement  of  T..  and  p ,  the  initial  (preshock)  test  gas  temperature 
and  pressure  respectively  and  U  the  ^incident  shock  velocity  is  in  principle 
most  straightforward  and  has  already  been  discussed  by  many  authors. 

However  in  the  present  experiments,  a  great  deal  of  care,  had  to  be 
taken  to  reduce  errors  in  the  measurement  of  these  parameters,  particularly 
as  they  define  the  boundary  conditions  for  calculating  the  shock  properties 
from  the  shock  Hugoniot  equations.  For  example,  a  1$  error  in  the  measurement 
of  U  generates  a  4$  error  in  the  calculated  value  of  N,  the  electron  density 
for  a  M  =20  shock  intp  1  mm  Hg  of  argon. 

S  ^  '  i  >.  • 

(i)  Initial  Preshock  Temperature,  T^.  >  .. 

The  initial  temperature  was  measured  by  means  of  a  mercury  buib  ther¬ 
mometer  calibrated  in  0.2°C  intervals  which  was  inserted  into  a  2  in.  deep  x 
l/4  in. dia.  oil  filled  port  located  close  to  the  test  section*  It  was  assumed 
that  the  test  gas  rapidly  comes  into  ,a  state  of  theimal  equilibrium  with  the 
shock  tube  in  the  time  period  between  test  gas  admission  and  shock  initiation 
(~  5  min).  Values  for  the  initial  temperature  were  interpolated  to  the  nearest 
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0.1°C  with  an  estimated  overall  error  of  +  0.1°C.  This  corresponds  to  a 
maximum  error  of  +  0.0%  in  the  absolute  value  of  the  initial  temperature. 

(ii)  Initial  (Preshock)  Pressure,  _ 

The  initial  pressure  p^  was  monitored  by  means  of  a;:  oil  manomete' 
based  on  the  design  of  Hayward  (Ref.  38) ,  which  covers  the  pressure  range  from 
0.1  to  bo  mmHg. 

The  manometer  consists  of  a  50  cm  long  x  i-cm  i.d.  U-tube,  which  is 
connected  at  one  end  to  an  oil  reservoir  and  at  the  other  end  to  the  vacuum 
system  via  a  length  of  flexible  vacuum  tubing.  By  rotating  the  manometer  through 
120°  about  a  horizontal  axis,  the  oil  (Dow  Corning  type  704  diffusion  pump 
oil)  is  allowed  to  collect  in  the  reservoir  and  is  vigorously  agitated  for 
several  minutes  so  as  to  induce  rapid  degassing.  Hayward  has  shown  that  this 
type  i> f  manometer  gives  a  much  better  correlation  with  a  McLeod  gauge-  than  the 
more  conventional  fixed  U  tube  oil  manometer  having  a-  cross  tube  with  stopcock. 

Dow  Corning  Corp.  (Ref.  39)  quote  a  specific  gravity  of  1.07  at  25°C 
for  DC704  oil  with  a  corresponding  temperature  coefficient  of  9*5  x  10“4/oq .  A 
laboratory  check  on  the  specific  gravity  of  the  sanple  used  produced  a  value 
of  I.069  at  25°C  which  was  therefore  employed  in  the  present  measurements  to¬ 
gether  with  the  temperature  coefficient  given  above. 

The  initial  preshock  pressure  p^  in  mmHg  is  obtained  from  the  difference 
in  oil  levels  H  mm  oil  by  the  relationship 
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where  p  o  =  13.595  is  the  specific  gravity  of  mercury  at  0°C  and  p  .  0  = 

1.069  sm  ^)  x  9*5  x  10  ^  is  the  specific  gravity  of  the  oil  at  T^SC  ’  ■L 
where  T^  is  the  laboratory  air  temperature  in  the  neighbourhood  of.  the  manometer. 

With  a  view  to  ensuring  an  accurate  measurement  of  p^,  the  manometer  was 
fitted  with  a  viewing  cursor  .and  back  view  mirror  so  as  to  prpvide  a  horizontal 
reference  mark  and  to  avoid  a  parallax  rading  error.  In  addition  the  vertical 
position  of  the  manometer  was  checked  from  time  to  time  with  a  plumb  line. 

Originally  it  was  only  intended  to  use  the  oil  manometer  for  pressures 
down  to  about  2  mmHg  and  to  use  a  McLeod  gauge  (5  nmHg  to  1  pHg)  for  lower 
pressures  down  to  6.3  mmHg.  However ,  in  the.  range  of  overlap  of  the  two  gauges, 
the  McLeod  gauge  gave  readings  which  were  about  %  higher  than  those  for  the 
oil  manometer.  This  was  initially  thought  to  be  due  to  ,a  poor  reference  vacuum 
in  the  isolated  limb  of  the  manometer  but  was  ruled  out  after  leak  testing  the 
manometer.  In  addition,  the  vapour  pressure  of  the  DC  704  oil  is  about 
2  x  10“ 8  mmHg  at  25°  C  (Ref.  39)  and  is  therefore  not  a,  contributing,  factor 
to  this  difference. 

In  the  absence  of  any  obvious  fault  in  the  oil  manometer,  the  error 
was  therefore  ascribed  to  the  calibration  of  the  McLeod,  gauge  (Stokes  Model 
No.  276-M) .  Consequently  all  p^  measurements  were  made  on  the  oil  manometer 

with  the  main  source  of  error  considered  to  be  that  due  to.  the  visual  reading. 
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of  the  oil  levels;  this  was  estimated  tobb.e  +  0.1  mm  of  oil  which  at  p^  =  0.3 
mmHg  corresponds  to  a  maximum  random  error  of  +  2§ $>. 

However,  in  order  to  attain  these  limits  of  accuracy"  at  low  p  values, 
the  oil  manometer  should  be  isolated  for  at  least  10  minutes  after  admission  of 
the  test  gas  before  H  is  read  in  order  to  allow  the  oil  adhering  to  the' tube 
walls  to  drain  back  through  the  meniscus.  A  ’settling'  reading  error  between  the 
left  and  right  hand  limbs  is  thereby  avoided. 

(iii)  Incident  Shock  Velocity,  U_  '  \  ’ 

* 1  1  ■  ■  ’  ’  ""  “  r ' '  1  s — 

The  incident  shock  velocity  is  obtained  by  the  usual  method  of  measuring 
the  time  interval  between  shock  arrival  at  two  axially  displaced  pressure  sensi¬ 
tive  detectors.  The  detectors  are  separated  by  609.6  +  O.i  mm  and. locatect  at 
points  which  are  equal  distances  upstream  and  downstream  of  the.  centre  of  the 
test  section  windows  (see  Fig.  2).  The  piezo-electric  shock  detectors  .(Atlantic 
Research  ID  25  Blast  Pressure  Transducers)  are  flush  mounted  and  have  6;  mm..dia  • 
lead  zirconate  titanate  elements  covered  with,  a  thin  translucent  layer  of  ^plasti¬ 
cized  epoxy  resin.  The  output  of  these  gauges  is  approximately  0.15.v/psl* 

The  signals  from  the  pressure  detectors  are  taken  directly  via  -short 
leads  to  battery-driven  emitter-follower  units,  having* unit  amplification,  and  ~ 
l/lOpsec  risetime  ard  then  to  the  start  and  "top  inputs  of  a  +  1/lOpsec  resblutiph 
electronic  timer  (Hewlett  Packard  Universal  counter  5325A).  This  counter- timer  , 
has  separate  variable  triggering  level  capability  for  each  channel. so  .that  varia¬ 
tions  in  sensitivity  between  gauges  <.an  be  accommodated.  However,  wherever- possible, 
matched  gauges  were  used  so  as  to  reduce  both  temporal  and  spatial,  errors;. incurred 
in  measuring  the  shock  velocity. 

The  finite  size  (6  mm  dia)  of  the  piezo-electric  elements  introduces1  a 
maximum  uncertainty  of  +  6  mm  into  the  distance  separating  the  gauges.  However, 
it  was  considered  that  by  using  matched  gauges,  this  uncertainty  could  be  reduced 
to  +  3  which  corresponds  to  +  l/2 $  uncertainty  in  the  gauge  separation.. 

The  risetime  of  the  output  waveform  of  these  gauges  in  response  to  the 
frozen  shock  front  pressure  jump  is  generally  of  the  order  of  lpsec  as  illustrated 
by  the  lower  trace  of  oscillogram  (b)  in  Fig.  3«-  The-  first. .  -ring'  of  these..  -, 
gauges  has  a  risetime  of  about  l/2  psec  and  an  amplitude  of  about  60$  of.  the  total 
pressure  rise.  The  separate  channels  of  the  timer  are  therefore,  set  to  trigger 
at  50$  of  the  anticipated  voltage  rise  for  the  frozen  shock  pressure  rise.  This 
introduces  a  maximum  uncertainty  of  +  1/2  psec  into  the  time  interval  measurement 
which  in  a  typical  period  of  100  psec  constitutes  a  random  error  of  +  l/2$. 

Hence  the  temporal  uncertainty  of  +  l/2 %  due  to  the  finite  rise  time 
of  the  gauge  combined  with  the  spatial  uncertainty  of  +  l/2 ijo  introduces  a  total 
possible  random  error  of  +  1%  into  the  measurement  of  the  shock  velocity.  Check 
measurements  made  using  two  pairs  of  similarly  displaced  gauges  with  a-,  similar 
counter  produced  discrepancies  between  the  two  readings  of  the  order  of  l/lO  psec 
indicating  that  the  +  1$  error  referred  to  above  is  somewhat  pessimistic. 

A  problem  encountered  with  these  pressure  gauges  at  an  early  stage  in 
the  project  is  illustrated  by  the  apparent  precursor  profile  in  -oscillogram-  (a) 
of  Fig.  3  which  tended  to  cause  the  various  electronic  units  to  pre- trigger.  - 

As  strong  shocks  into  argon  are  known  to  produce  significant  photon  and 
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electron  precursors  (Ref.  4o),  it  was  suspected  that  the  precursor  signal  monitored 
by  the  pressure  gauges  was  due  to  one  or  both  of  these  aforementioned  sources. 


In  the  first  instance  the  translucent  epoxy  resin  covering  the  gauge  was 
given  an  additional  covering  of  transparent  silicpne  (General  Electric  RTV)  rubber. 
This  had  no  effect  in  blocking  the  precursor  signal  suggesting,  that  electron  leak¬ 
age  to  the  gauge  elements  was  not  the  cause.  However,  a,  similar  film,  of  black  sili¬ 
cone  rubber  completely  blacked  this  precursor  as,  illustrated  by  oscillogram  (b) 
in  Fig.  3j  suggesting  that  the  precursor  was  due  solely  to  the  photon  flux  coming 
from  the  strongly  radiating  equilibrium  region  behind  the  shock  front.  This  is 
not  surprising  as  it  is  well  known  that  the  piezo-electric  effect  is  always,  accom¬ 
panied  by  the  pyro-electric  effect  in  which  an  electric  charge  is  produced  at 
crystal  surfaces  exposed  to  thermal  radiation  (Ref.  31) • 

The  general  consensus  in  the  literature  concerning  the  precursor  phenomenon 
is  that  for  strong  shocks  intp  inert  gases  the  observable  precursor  electron  den¬ 
sity  is  produced  by  photoionization  (Ref.  4(5)  whereas  for  relatively  weak,  shocks, 
e.g.,  Mg<  10  into  argon,  the  observed  precursor  electrons  are  now  thought  to  be 
due  to  electron  diffusion  (see  Refs.  40,  42,  43)  which  for  conditions  of  high  Mach 
number  becomes  masked  by, the  photoionization  process. 

It  is  felt  that  some  attention  should  be  paid  to  the  effect  of  shock 
attenuation  both  on  the  shock  velocity  measurement  and  on  the  plasma  properties 
behind  the  shock  front.. 

For  shock  Mach  numbers  from  17  to  24  and  initial  pressures  from.  0.3  to 
3  mm  Hg,  the  shock  front  attenuation  in  the  .region,  .of  the  test  .section  was  very 
low  being  typically  1/3$  per  foot  of  the  7  in.  x  4  in1,  shock  tube..  The,  attenua¬ 
tion  as  monitored  by  several  velocity  measurements  in  the  test  section  region  was 
generally  so  low  as  to  be  less  than  the  uncertainty  in  the  measurement  of  the 
shock  velocity  itself  which  has-been  shown  to  be  .of  the  order  +  1 

As  will  be  indicated  in  Section  4,  the  ionization  relaxation,  region  from 
the  shock  front  to  conditions  of  Saha  equilibrium  is  usually  of  -the  order  of 
10  to  20  mm  so  that  the  effect  of  attenuation  on  the  shock  properties  over. this 
distance  can  be  neglected  in  view  of  the  attenuation  .figures  given  above. 

3.2  Optical  Diagnostics  and-  Recording  Equipment. 

The  description  of  the  optical  system  is  broadly  divided  into  arc^s  deal¬ 
ing  with  the  laser  light  source,  the  interferometer  itself  and  finally  the  photor 
graphic  recording  system. 

3.2.1  Laser  Light  Source 

Dual-frequency  lasers  have  been^  employed  by  a  number  of  authors,  to- 
observe  various  plasma  phenomena  either  by  interferometric  or  schlieren  .techniques 
(Ref.  44  to  49). 

The  dual  frequency  laser  employed  is  a  TRG  104a  system,  which  is  giant 
pulsed  by  means  of  an  integraliy-mounted  Pockels  cell  Q- switch  and  is  capable  of 
operating  either  as  a  ruby  or  a  neodymium  doped  glass  laser  by  suitable  choice 
lasing  element  and  cavity  end  reflectors.  The  ruby  laser  operates  at  a  wavelength 
of  6943. 5A  at  24 °C  (Ref.  80)  with  a  spectral  bandwidth  of  less  than  O.lAand  when 
Q-switched  produces  a  30-Mtf  pulse  of  15-nsec  duration.  Similarly  the  neodymium 
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laser  operates  at  IO600X  with  a  bandwidth  of  about  50$  for  which  the  Q-switched 
power  and  pulse  duration  are  20-MW  and  15-nsec,  respectively. 

The  second  harmonic  generator  for  the  ruby  laser  consists  of  a  critically 
oriented  crystal  of  potassium  di-hvdrogen  phosphate  (KDP)  which  is  mounted  in  an 
oil-filled  cell  attached  directly  ■  the  output  end  of  the  laser  cavity.  Power 
generation  for  this  process  is  typically  5  to  10$  efficient  for  the  3,-cm-  long 
crystal  and  produces  coherent  radiation  at  3%1.75$>  which  is  both  parallel  to 
and  concurrent  with  the  incident  laser  radiation.  The  spectral  bandwidth  of  the. 
second  harmonic  of  the  ruby  laser  is  typically  0.05 A  or  less  .  Similarly  -an 
appropriately  cut  KDP  crystal  is  used  to  achieve  second  harmonic  generation,  of 
the  neodymium  laser  output  at  5300$  with  a  bandwidth  of  less  than-  25$, 

A  certain  measure  of  care  should  be  taken  in  the  adjustment  and  use  of 
these  second  harmonic  generators.  Although  the  fundamental  and  second  harmonic 
beams  are  parallel,  the  second  harmonic  beam  becomes  displaced,  from. ‘the  axis*  .of 
the  fundamental  due  to  double  refraction  within  the  KDP  (Ref.  50).  For  both  the 
ruby  and  neodymium  lasers,  the  angular  displacement,  of  the  second  harmonic  in 
KDP  is  approximately  1.77°  (Refs.  50,  51) >  which  for  the  3-cm  path-length  crystal 
currently  employed,  corresponds  to  an  output,  beam  displacement  of  0.92  man  It  is 
therefore  clear  that  when  the  dual  frequency  is  used-  to  illuminate  a  schlieren 
system,  the  laser  should  be  arranged  so  that  the  axes  of  both  the  fundamental 
and  second  harmonic  beams  intersect  the  schlieren,  knife  edge, 

The  efficiency  of  second  harmonic  generation  is  critically  dependent 
on  the  orientation  of  the  KDP  -crystal  with  respect  to*  the  incident  laser  beam 
(Ref.  52)  necessitating  careful  adjustment.  Second  harmonic  generation  is  also 
strongly  temperature,  dependent  (Refs.  52,  53)  such  that- small  changes  in  the 
tempe  rature  of  the  KDP  induced  by  variable  laboratory  conditions  or  by  frequent 
laser  operation  can  result  in  a  large  reduction  in  tsecond-harmonic  power  .conver-r 
sion.  Ideally  this  unit  should  be  contained  in  a  thermostatically  controlled 
environment.  For  second-harmonic  .generation  by  the  ruby  laser  (Ref,  53)  and 
the  neodymium  laser  (Ref.  52)  in. KDP,  the  temperature  gradients  for  the  phase, 
matching  angle  are  46arfisec/°C  and.  33.6  arcsec/C  respectively. 

In  the  field,  of  optical  plasma  diagnostics,,  the  advantages  of  -the 
pulsed  laser  over  the  more  conventional,  spark  or  exploding  .wire,,  light  are  -as . 
follows:  ,  ' . 

(i)  The  short  pulse  duration  of  the  Q-switched  laser,  generally  from  -10  to 

30  nsec,  is  about  an  order  of  magnitude  shorter  than  for  thermal  sources, 
of  adequate,  intensity..  For  the  laser.,  this  results  in  :exceilent  shock 
front  and  relaxation  zone  definition  when  photographic  recording,  is» 
employed.  ,  .  .  , 

(ii)  The  high  intensity  laser  radiation  is  ideal  for  discriminating  against 
plasma  radiation  particularly  when  employed  in  conjunction,  with-.a  com¬ 
bination  of  spatial  and  spectral  filters.  However,  care.  must,  .be  .takeh, 
to  (insure  that  the  spatial  power  density  of  the  laser,  -radiation  is 

not  sufficient  to  cause  breakdown  of  the  transmission  or  reflection 
media.  If  the  beam  is  brougb ;  to  a  sharp  focus.. either  in  air  or  some- 
other  gaseous  system,  local  plasma  formation  “ay  occur  resulting  in 
a  loss  of  usable  laser  power  due  to  absprption.-  In^the  case  of  mirror 
and  filter  surfaces  or  transmission  media  such  as  absorption  filters,  . 
lenses,  windows,  etc.,  permanent  radiation  damage  can  occur. 
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(iii)  The  high  degree  of  spectral  (temporal)  coherence  of  the  laser  output 
makes  for  excellent  wavelength  definition  e.g.}  the  ruby  laser  wave¬ 
length  at  2^°C  is  69^3. 5$  (Kef.  80)  with  a  bandwidth  of  0.1$  when 
resonant  cavity  end  reflectors  are  used.  With  a  white  light  source, 
an  acceptable  degree  of  monochromaticity  can  be  obtained  by  means  of 
interference  filters  but  at  the  expense  of  a  large  drop'  -in  intensify. 

As  discussed  in  Section  2.1,  it  is  necessary  for  the  laser  wave¬ 
length  to  be  remote  from  any  absorption  lines  or  bands  of  the  trans¬ 
mitting  medium  under  investigation,  in  order  to  avoid  the  effects  of. 
anomalous  dispersion  or  absorption.  However,  in  certain  circumstances, 
these  effects  can  bei'utilized  to  selectively  monitor  the  populations 
and  behaviour  of  specific  internal  states ,  particularly  those  involved 
in  transitions  which  control  a  given  relaxation  or  reaction  process. 

(iv)  The  ability  of  the  pulsed  laser  to  operate  over  a  wide  range  of 
wavelengths  -from  the  near  ultra-violet  (3^71$)  to  the  near  infra¬ 
red  (10600A)  makes  it  an  ideal  light  source  for  plasma  interfero¬ 
metry.  This  laser  therefore  operates  over  the  full  spectral  band¬ 
width  of  the  present  optical  system.  The  second  harmonic  of  the 
the  ruby  laser  at  3^71$  lies  just  below  the  cutoff  frequency  of 
conventional  glass  optics  whereas  the  fundamental  of  the  neodymium 
laser  at  10600A  lies  -at  the  long  wavelength  sensitivity  limit  of 
currently  available  photographic  emulsions. 

The  two-wavelength  plasma  interferometry  techniques  of  -Alpher 
and  White  (Ref.  2)  for  measuring  electron  density^  requires  a  long 
wavelength  light  source  operating  simultaneously-  at  two  different 
wavelengths  (see  Section-  2 . 3) i  This  role  is  ideally' met  by  the 
neodymium  laser  and  its  second1  ’harmonic  generator-  operating  simul¬ 
taneously  at  10600$  -and  5300$-s 

For  the  observation  of  changes-  in'  mass  density  a  short  wave¬ 
length  source  is  required  (see  Section-  2.2)  and- is  accomplished  by 
use  of  the  ruby  laser  second  harmonic  at  3^7'1A. 

(v)  By  the  nature  of  its  formation,  the-  laser  beam  is  highly  unidirectional, 
being  in  the  form  of  a  narrow  hear  parallel  beam  with  a  divergence  of 

a  few  milliradians  such  that  a  system  of  condensing  optics  is  not 
necessary. 

- 

(vi(t»  Due  to  the  monochromatic  and  unidirectional  state-  of  the  laser  beam', 
it  can  (in  principle)  be  focussed  down  -to  a  diffraction  limited  spot. 

The  laser  is-  therefore  ah  ideal  source  for  use  in  schlieren  photo-, 
graphy  as  the  sensitivity  of  this  technique  is  an  inverse  function 
of  source  size  (Ref,  82). 

However,  plasma  formation  and  knife  edge  damage  may  occur'  place- ' 
the  laser  power  exceeds  some  critical  value-.  Al-so,  fine  control  of 
the  small  source  image  on  the  knife  edge  requires  ah.  extremely  fine- 
degree  of  control. 

(vii)  The  high  degree  of  spatial  and  temporal  'coherence  of  the  laser  beam 
makes  it  an  ideal  interferometer  light  source.  When  a  white  light 
(thermal)  source  is  used  to  illuminate  a  Mach-Zehnder  type  interferometer. 
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it  is  necessary  to  path  length  match  the  two  beam  paths  to  within  half 
a  wave  (A/2)  if  a  fringe  field  of  high  contrast  is  tq:.be  obtained. 
However,  with  the  laser,  this  requirement  can  be  relaxed. 

If  the  bandwidth  of  the  laser  output  is  AA  for  laser  radiation, 
of  wavelength  A,  tnen  it  can  be  easily  shown  by  a  classical  wave 
optics  argument  that  the  allowed  path  length,  mismatch  over  which  the 
two  beams  become  A/2  out  of  phase  is  given  by  A^/2AA<j  For  the  r.uby 
laser  operating  at  6943$  with  a  spectral  bandwidth  of  0.1$,  this  means 
that  the  path  length  of  the  two  beams  must  not  differ  by  more,  than 
2§  cms  if  high  fringe  contrast  is  to  be  obtained,  --The^e odymium  laser 
operating  at  10600A  has  a  bandwidth  of  about  50$  so' that  in  this 
case  a  path  length  mismatch  of  only  0.1  mm  is  allowed. 


However,  the  direct  use  of  the  laser  beam  .does  carry  certain- dis¬ 
advantages.  Firstly,  the  spatial  or  transverse  modes  of  the.  laser  beam 
result  in  non-uniform  illumination  of  the  syst  jn  under  investigation 
unless  spatial  mode  selection  has  ..been  employed  usually  resulting  in  a 
reduction  in  beam  intensity.  Secondly,  due  to  the  coherent  .nature-  of 
the  laser  radiation,  small  imperfections  and  dust  particles  in  the 
optical  system  tend  to  generate  a  system  of  diffraction  rings  which 
are  clearly  visible  in  the  photographic  plane,  finally;  numerous 
sets  of  unwanted  fringes  tend  to  form  by  interference  between,  the  many 
optical  surfaces  in  the  system  all  of  which  become  visible  in  -the 
photographic  plane.  As  will  be  described  in  Section  3  •, 4,  these  diffir 
culties  can  be  eliminated  by  spoiling  the  phase  coherence  of  the  laser  . 
beam  by  transmission  through  a  diffuser,  which  in  effect  becomes  the 
interferometer  light  source. 


3.2.2  Electronic  Operation 


The  operation  of  the  laser  as  the.  shock  tube=?interferometer  iight  source, 
is  complicated  by  the  fact  that  the  lasing,  elements  must  be  flash  lamp  pumped 
for  a  fixed  period  of  time  prior  to  Q-switching  in  order  that  the  required, 
population  inversion  is  attained. 


The  Pockels  cell  Q-switch  is  arranged  to  open  a,t  a  -time  when  the  flash 
lamp  output  is  just  beginning  to  decay.  If .  the  Pockels,  cell  is-  opened  at,  an. 
earlier  time,  then  in  addition  to  the  required  giant  pulse,  a  series  of  normal 
mode  pulses  will  be  obtained,  as  the  quarterr-wave  voltage  is  re-applied  ,to  the 
Pockels  cell  over  a  period  of  about  1  msec..  This  will  result  in  a.  number  of 
undesired  interferograms  being  super-imposed  on  the  desired  interferogram  of 
the  shock  phenomenon  under  investigation.  Alternatively,,  if  the.  Pockels  -cell 
is  opened  at  a  later  time,  then  a.. much  weaker  giant  pulse  will  be  obtained 
(nef .  54) .  The  flash  lamp  pumping  times  for  the,  ruby  and  neodymium  lasing  ele-: 
meats  as  used  in  the  TEG  104A  laser  system  .are  950  +  25 -usee  and  -80P  +  20  psec 
respectively.  The  800  usee  delay  period  for  the  neodymium  laser  provides  a 
giant  pulse  which  is  just  sufficient  for  interference  fringe  photography  but 
which  subsequently  allows  the  laser  to  free  run  in  normal  mode  operation.  This 
effect  manifests  itseif  as  a  loss  in  fringe  resolution  rather  than  .as  a  series 
of  multiple  exposures  as  can  be  seen  in  the  interferogram  produced  by  the  neo¬ 
dymium  laser  shown  in  Fig.  l?(ii,).  Howeyer  the  second  har  mimic  generation  process 
varies  as  the  square  of  the  fundamental  laser  pulse  power  so.  that  these  normal 
mode  pulses  are  too  weak  to  generate  a  perceptible  level  of  second  harmonic 
radiation.  This  accounts  for  the  relatively  'noise*  free  interferogram  produced 


by  the  second  harmonic  shown  in  Fig.  15(i).  Fortunately  the  ruby  laser  can  be 
Q- switched  after  a  longer  delay  period  (  ~  950  psec)  without  incurring  a  serious 
loss  in  giant  pulse  power  but  where  the  flash  radiation  has  decayed  to  a  level 
which  is  below  threshold  for  normal  mode  laser  operation. 

Hence  a  certain  knowledge  of  the  shock  velpcity  and  attenuation  is 
required  prior  to  each  shock  run  in  order  that  the  laser  flash  iamp  can  be  arranged 
to  trigger  at  the  appropriate  time  prior  to  shock  arrival  in  the  test  section. 

A  schematic  of  the  system  by  which  this  is  achieved  is  shown  in  Fig.  2.  Shock 
arrival  at  pressure  gauge  1  triggers  the  delay  generator  which  after  the  elapse 
of  the  delay  period  triggers  the  flash  lamp  capacitor  bank.  The  delay-generator 
period  plus  the  flash  lamp  pump  period  is  therefore  equal  to  the  time  the  shock 
front  takes  to  travel  from  detector  1  to  the  centre  of  the '.test Lsection. 

The  Pockels  cell  Q,-sv/itch  is  similarly  activated  via  a  delay  circuit 
initially  triggered  by  shock  detector  2.  The  delay  period  in  this  case  is 
arranged  to  be  such  that  the  Pockels  cell  opens  when  the  shock  front  is  in  the 
desired  position  in  the  test  section  aperture.  The  delay  generator  for  the 
Pockels  cell  is  conveniently  triggered  by  the  synchronized  output  signal  of  the 
'start'  channel  of  the  counter  timer. 

3.2.3  Mach-Zehnder  Interferometer 

Tile  construction  and  operation  of  the  9  in.  dia.  aperture  Mach-Zehnder 
interferometer  has  been  described  by  Hall  (Ref.  55).  The  adjustment  procedure 
for  obtaining  and  focussing  the  fringes,  which  is  ^essentially  the  same  as  described 
by  Hall  but  with  small  refinements  and  modifications  resulting  from  the  installation 
of  the  pulsed  laser  light  source,  will  be  described  in  the  following  section 
(3.2.4). 


A  scale  drawing  of  the  interferometer  is  shown  in  -Fig.  4.  The. .central 
body  of  the  interferometer  is  formed  by  a  3/8  in.  thick  welded  U-shaped  box 
frame  structure.  The  projections  at  either  end  supporting  the  parabolic  colli¬ 
mating  and  condensing  mirrors  are  canvas  coyered  to  prevent  dust  and  air  currents 
from  entering  the  interferometer. 

During  operation  of  the  interferometer,  it  was  found  that  small  changes 
in  the  laboratory  temperature  had  a  marked  effect  on  the  fringe  settings.  The 
small  rise  in  the  laboratory  temperature  generated  by  the  electronic  equipment- 
used  in  the  shock  tube  experiments  -was  found  to  have  a  considerable  effect  on  the. 
setting  of  the  focus,  orientation  and  spacing  of  the  fringes  during  the  4  to 
5  minute  interval  between  the  taking  of  the  no-flow  and  flow  interferograms.. 

Further  investigations  shoWed  that  the  steel  shell  of  the  interfero¬ 
meter  was  expanding  or  contracting  almost  instantaneously  in  response  to  .these 
small  changes  in  the  laboratory  air  temperature.  Temperature  changes  of  the 
order  of  l/lOO°C  were  found  to  induce  changes  in  the  interferometer  setting 
sufficient  to  generate  movements  of  the  order  of  l/lO  fringe  at  3471&. 

By  lagging  the  outside  of  the  interferometer  shell  with  1  in.  thick 
foam  sponge,  the  response  time  of  the  interferometer  to  temperature  changes 
was  increased  from  about  30  seconds  to  about  TO  minutes  Thereby  providing  an 
adequate  delay  period  in  which  to  conduct  the  shock  tube- interferometer  experi¬ 
ments.  Additional  precautions  involved  placing  non- interferometric  quality 
windows  over  the  entry  and  exit  parts  of  the  interferometer-  and  enclosing' the 
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optical  paths  between  the  shock  tube  test  section  and  interferometer  with  foam 
sponge  tubes  with  a  view  to  preventing  the  entry  of  both  dust  particles  and  air 
currents  into  the  interferometer. 


The  laser  system,  the  interferometer,  the  test  and.  compensating  chamber 
sections  and  the  camera  are  all  initially  centered  and  co-aligned  by  means  of  a 
heli,um-neon  laser  as  shown  in  Fig.  5(a).  The  purpose  of  the  pentagonal  prism 
is  to  ensure  that  the  optical  axis  of  the  helium-neon  laser  effectively  coincides 
with  that  of  the  ruby/neodymium  laser. 

As  the  position  of  the  test  section  windows  is  fixed.,  all  other  components 
in  the  system  including  the  compensating  chamber  must  be  positioned  with  respect 
to  the  test  section  aperture.  The  optical  axis  of  the  interferometer  is  arranged 
to  be  perpendicular  to  the  shock  tube  axis  by  adjusting  the  various  system  com¬ 
ponents  until  the  reflection  of  the  laser  beam  from  the  test  section  windows  is 
returned  to  the  laser  cavity. 

As  discussed  in  Section  3.2.1,  accurate  path  length  matching  of  a  two 
beam  interferometer  is  not  necessary  when  using  a  laser  light  source  having  a 
long  coherence  length.  However,  this  leads  to  the.  undesirable  effects  of  laser 
mode  structure,  spurious  fringes  and  dust  generated  diffraction  patterns  being 
visible  in  the  photographic  plane.  These  effects  can  be  eliminated  by  passing 
the  laser  beam  through  a  diffuser  which  becomes  in  effect  a. light  source  of  finite 
size.  The  spatial  modes  of  the  laser  beam  become  smeared. out  resulting  in  uniform 
illumination  of  the  test  section  whereas  the  diffraction  rings  and  the  fringes  are 
effectively  eliminated  due  to  the  reduced  spatial  coherence  of  the  beam  transmitted 
by  the  diffuser.  The  diffuser  is  therefore  located  at  the  focus  of  the  collimating 
optics  of  the  interferometer,  which  in  this  case  is  a  60  in.,  focal  length  para¬ 
bolic  mirror  (see  Fig.  4).  As  a  result  of  this  expedient,  it  now  becomes  neces¬ 
sary  to  path  length  match  the  interferometer  according  to  the  procedure  described 
by  Hall  (Ref.  55).  ’ 

Path  length  matching  is  initially  achieved  using  a  high  pressure  mercury-? 
vapour  lamp  light  source  together  with  a- condensing  lens  and- source  slit  as  -shown 
in  Fig.  5(b)  (see  Ref.  55  Tor  details  concerning,  adj  ustment  of  interferometer 
optics*').  Once  fringes  are  visible  on  the  camera  screen,  the  central  order  fringes 
can  be  located  using  a  hand  spectroscope.  By  imaging  three  or  four  horizontal 
fringes  onto  the  vertical  spectroscope  slit,  oblique  fringes  are  observed  in 
the  spectroscope.  Beam  splitter  1  (see  Fig.  4)  is  then  translated  until  the 
fringes  in  the  spectroscope  become  horizontal.  These  are  then  the  central 
order  fringes  where  the  two  optical  paths  of  the  interferometer  are  equal  for 


*  Hall  (Ref.  55)  omitted  to  mention  in  his  procedure  that  after  path  length 
matching  of  the  interferometer  by  translation  of  beam  splitter  1  as  described 
in  Section  3.2  of  Ref.  55?  the  two  beams  are  ho  longer  exactly  superimposed. 
This  therefore  necessitates  further  adjustments  of  mirrors  1  and  2  as  des¬ 
cribed  in  Section  3.1  of  Ref.  55  although  this  iterative,  procedure,  rapidly 
converges  to  point  where  high  contrast  monochromatic  fringes  can  be  seen. 


all  wavelengths  to  within  a  distance  of  the  order  of  one  wavelength. 

Both  the  virtual  fringes  focussed  in  the  test  section  and  the  image  of 
the  test  section  are  then  brought  to  a  real  focus  on  the  camera  screen  by  use 
of  a  large  diffuse  light  source  (see  Ref.  55  for  procedural  details).  This  con¬ 
veniently  takes  the  form  of  a  domestic  tungsten  lamp  as  shown  in  Fig.  5  (d). 

The  pulsed  (ruby  or  neodymium)  laser  is  then  centered  and  aligned  with 
the  optical  axes  of  the  system  using  a  helium-neon  laser  and  pentagonal  prism 
as  shown  in  Fig.  5(e)* 

The  optical  system  as  employed  in  the  dual  frequency  laser  interfero¬ 
meter  experiments  is  shown  in  Fig.  5(f).  The  1-cm  dia.  dual-frequency  laser 
beam  is  reduced  to  a  spot  of  about  5~vm  dia.  bn  a  double  sided  diffuser  by  means 
of  a  weak  single  component  lens.  The  diffuser  was  produced  by  grinding,  both 
faces  of  a  microscope  slide  with  a  10  p  alumina  or  carborundum  grit.  The .effective 
source  size  was  then  defined  by  a  5-nm  dia.  stop  affixed  to  the  output  side  of 
the  diffuser  (see  Fig.  5(f)). 

Although  use  of  this  diffuser  resulted  in  a  considerable  reduction  in 
the  intensity  of  the  interferometer  beam,  it  was  still  found  to  be  sufficient 
to  activate  300  ASA  panchromatic  emulsions  at  69I3&  while  simultaneously  pro¬ 
viding  adequate  discrimination  against  the  argon  shock  plasma  radiation. 

In  principle  the  finite  size  of  the  source  aperture  tends  to  degrade 
the  degree  of  cbliimation  of  the  interferometer  beam  and  also  to  limit  the  .number 
of  fringes  visible  in  the  plane  of  focus.  For  the  present  system,  the  5-mm 
diameter  aperture  corresponds  to  an  interferometer  beam  divergence  of  about 
0.1°  which  over  the  1.010-in.  width  of  the  shock  tube  path  length,  results  in  a 
limiting  spatial  resolution  of  about  0.01  in.  Similarly  the  number  of  visible 
fringes,  given  by  l/02  (Ref.  56).  where  0  is  the  interferometer  divergence  half 
angle,  is  of  the  order  of  105.  It  is  therefore  clear  that  the  finite  size,  of 
the  diffuser  light  source  does  not  constitute  a  real  limitation  on  the  use  of 
this  system  in  the  present  circumstances. 

3.2,5  Test  Section  and  Compensating  Chamber  -Windows 

The  8  in.  dia.  aperture  x  3-3/1  in.  uhick  test  section  and  compensating 
chamber  windows  were  manufactured  from  Schott  BK7  optical  quality  glass  and 
were  polished  flat  so  that  the  8  in.  diameter  transmitted  wavefront  of  each 
window  was  flat  to  A/1  at  6328A  with  wedge  angles  not  exceeding  1  sec  of  arc. 

The  windows  were  mounted  so  as  to  provide  a  disturbance  free  flow 
toQthe  shock  wave.  However,  during  manufacture,  the  inner  edge  was  given  a 
15  x  1/8"  wide  bevel  as  a  precaution  against  chipping.  Although  this  de¬ 
pression  did  not  generate  any  interferometrically  observable  flow  disturbances, 
the  flow  stagnation  in  this  region  had  the  effect  of  causing  severe  erosion 
giving  the  glass  surface  a  sand  blasted  appearance.  This  effect  was  satisfac¬ 
torily  eliminated  by  filling  the  bevel  in  with  flexible  type  of  epoxy  resin. 

In  order  to  facilitate  removal  of  the  windows  from  the  frames,  the  bevel  surface 
was  treated  with  an  anti-bonding  agent  so  that  the  epoxy  resin  only  adhered  to 
the  window  frames. 

A  serious  problem  with  these  test  section  windows  as  with  a  previous 
pair  of  rectangular  aperture  windows  (Ref.  35)  involved  the  appearance  of  small 
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cracks  originating  from  the  re-entrant  corners  of  the  *T*  shaped  cross  section  of 
the  windows.  These  cracks  were  found  to  he  caused  by  the  application  of  uneven 
clamping  pressure  rather  than  by  shock  wave  pressure  loading.  Although  etching  of 
the  re-entrant  corner  region  with  hydrofluoric  acid  is  known  to  remove  surface 
microcracks  from  which  these  cracks  originate,  optical  laboratories  were  not  pre^- 
pared  to  perform  this  work  in  view  of  the  proximity  of  the  interferometric  quality 
surfaces.  Future  window  designs  have  therefore  abandoned  the  •T*  shaped  cross 
section  support  mounting  for  a  bevel  edged  support  format; 

3.2.6  Interference  Fringe  Photography 

The  camera  used  for  performing  two  wavelength  interferometry  is  shown  in 
plan  view  in  Fig.  it  and  schematically  in  Fig.  5(f). 

The  interferometer  beam  is  brought  to  a  point  focus  at  the  iris  diaphragm 
which  acts  as  a  stop  for  the  uncollimated  plasma  radiation.  The  image  of  the  test 
section  is  then  focussed  onto  the  photographic  plane  by  an  80  cm  focal  length  con¬ 
verging  meniscus  type  spectacle  lens  located  just  behind  the  iris  diaphragm. 

An  aluminized  beam  splitter  was  found  to  be  adequate  for  separating  the 
two  beams  although  in  circumstances  of  low  intensity  the  use  of  a  colour  selective 
dichroic  beam  splitter  would  be  more  advantageous.  Narrow  band  interference  filters 
centered  on  the  respective  laser  wavelengths  provide  further  discrimination  against 
the  plasma  radiation.  It  is  worth  noting  that  in  none  of  the  experiments  to  be 
described  in  Section  4,  was  the  plasma  radiation  sufficiently  intense  as  to  be 
observed  on  the  interferograms .  This  contrasts  sharply  with  the  experience  of 
Alpher  and  White  (Ref.  2)  who  employed  a  spark  light  source  to  perform  two  wave¬ 
length  interferometry  on  shock  generated  argon  plasmas. 

Due  to  the  dispersive  nature  of  the  camera  lens,  care  had  to  be  taken  to 
ensure  that  the  photographic  emulsion  was  correctly  positioned  for  each  of  the 
four  laser  wavelengths  employed  In  these  experiments .  Precise  focussing  of  the 
centre  of  the  test  section  in  the  photographic  plane  was  accomplished  by  mounting 
a  grid  of  0.010  in.  diameter  wires  on  both  sides  of  the  test  section  windows-.  A 
sharp  focus  of  the  .test  section  was  then  obtained  by  moving  the  camera  viewing 
screen  to  a  point  midway  between  the  position  of  focus  for  the  two  sets  of  wires. 

As  these  wires  were  also  used  as  reference  markers  from- which  the  fringe. shifts 
were  measured,  care  was  taken  to  ensure  that 'the  wires  lay  either  parallel  to  or 
perpendicular  to  the  shock,  tube  axis.  For  the  present  optical  system,  the  depth 
of  field  was  sufficiently  large  that  both  sets  of  wires  remained  in  reasonably 
sharp  focus  when  the  camera  was  focussed  onto  the  centre  of  the  test  section. 

Eastman  Kodak  RS  Panchromatic  (650  A.S.-A)  emulsion  on-  0.007  in.  thick 
estar  sheet  base  was  found  to  provide  good  response  to  the  laser  radiation  at 
3471  X,  5300  2  and  6943  X  when  processed  in  Kodak  D19  developer. 

* 

Interference  fringe  photography  at  10600  %,  the  wavelength  of  the 
neodymium  laser,  presents  more  of  a  problem.  Eastman  Kodak  I-Z  spectroscopic, 
plates  for  use  in  the  near  infrared  region  were  initially  used  in  the  unsensitized 
condition  but  failed  to  give  a  sufficiently  dense  image .  However  when  sensitized 
in  chilled  ammonia  solution  immediately  prior  to  use  according,  to  the  prescription 
of  Pope  and  Kirby  (Ref.  57),  the  speed  of  the  I-Z  emulsion  is  increased  by  a 
factor  of  about  400  which  was  more  than  adequate  for  the  purposes  of  the  present 
experiments . 
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A  final  check  on  the  interference  fringe  field  setting  for  the  visible  • 
and  ultraviolet  laser  wavelengths  was. made  immediately  prior  to  a  shock  run  using 
Polaroid  (3000  A.S.A.)  positive  pript  film  thereby  ensuring  the  satisfactory  funr- 
ctioning  of  the  optical  -  electronic  system. 

3.3  Data  Analysis 

The  procedure  adopted  for  extracting  the  fringe  shift  data  from  the  flow 
and  no-flow  interferograms  is  described  in  Appendix  A.  However  there  remains  the 
problem  of  making  a  judgement  with  x-egard  to  the  location  of, the  centre  of  any 
given  fringe. 

Preliminary  tests  with  a  microdensitometer  showed  that  the  dark  fringes 
gave  relatively  sharp  extrema  which  were  easy  to  locate.  Due  to  the  over-exposure 
of  the  bright  fringes,  these  fringes  become  considerably  broadened  such  that. the 
corresponding  dark  fringes  in  the  limit  of  very  high  intensity,  tend  to  become 
well  defined  thin  lines. 

However  a  sufficiently  accurate  location  of  the  dark  fringe: scentres 
was  obtained  by  enlarging  the  interferograms  by  a  factor  of  10  and  making  a  direct 
visual  judgement  of  fringe  location.  The  enlargements  were. made  on  Kodak  1594 
Aerial  Mapping  Paper  which  maintains  a  high  degree  of , dimensional  stability  after 
processing.  It  also  has  a.  matt  (i.e.,  non-glossy)  finish  which  makes  for  easy- 
visual  interpretation. 

The  location  of  the  dark  fringe  centre  lines  from  the  positive  enlarge¬ 
ments  were  made  to  an. estimated  accuracy  of  less-  than  +.0.2  mm.  This  corresponds 
to  a  maximum  error  of  +  0 . 04  fringes-  at  34.71.  ft  and  5300  X  for  a  minimum,  fringe 
spacing  cf  10  mm  and  +  0.02  fringes  at  6943  ft  and  10600  ft  for  a  minimum  fringe 
spacing  of  20  mm. 

4.  EXPERIMENTAL  MEASUREMENTS ,  RESULTS,  DISCUSSION  AND  CONCLUSIONS 
4.1  Introduction 

,  , - -  ■  ,r  -—I . 

The  structure _of  the  relaxation,  equilibrium  and  recombination  -regions 
behind  strong  shocks  into  argon  has  already  been  examined  both  experimentally  and 
theoretically  by  Bershader  and  his  c.o-workers-.  (Refs..  58,  59  >  60).  Hence  onjy  a 
brief  qualitative  description  of  the  non-equilibrium  processes  is  given  here  in 
order  to  familiarize  the  reader  with  the  context  t>f  the  -present .  experiments-; 

Figure  6(i)  shows  typical  fringe  shift  profiles  behind  a  strong  shock 
into  argon  for  both  the  laser  fundament  al(?\)  and  second  harmonic  (A/2)  wavelengths. 
The  positive  fringe  shifts  across  the  shock  front,  is  due  to  the  compression  of, 
the  neutral  atom  where  the  shift  for  the  laser  second  harmonic  is  approximately 
twice  that  for  the  fundamental  wavelength  due  to  the  inverse  linear  dependence 
of  the  heavy  panicle  fringe  shift  on  wavelength  (see  Sqe t?2223).  On  approaching 
equilibrium  there  is  a  rapid  increase  in  the  free  electron  concentration,  suffi-  : 
cient  to  produce  a  large  negative  fringe  shift  particularly  at  the  laser  funds- 
mental  wavelength;  this  is  due  to  the,  linear  dependence  of  the  free  electron 
fringe  shift  on  wavelength  (see  Eq.  2.1.5-  f-nd  Eq.  2.2.3). 

Figure  6(ii)  demonstrates  the  two  temperature  nature  of  the  incubation 
zone  prior  to  equilibrium  where  the  neutral  atoms  are  initially  brought  to  a 
relatively  high  translational  or  'frozen'  temperature  immediately  behind  the 
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shock  front.  This  gas  of  highly  energetic  neutral  atoms  then  acts  as  an  energy 
source  both  for  the  ionization  process  and  for  the  thermalization  of  the  free 
electrons . 

Finally  Fig.  6(iii)  illustrates  in  terms  of  the  changes  in  the 
degree  of  ionization  the  various  stages  of  non- equilibrium  which  exist  behind  a 
strong  shock  moving  into  an  inert  gas.-  Regime  I  (in  the  nomenclature  of  Ref. 

58)  is  the  region  in  which  ionization  occurs  through  the  relatively  inefficient 
atom-atom  collision  process  in  the  absence  of  a  significant  concentration  of  free 
electrons.  Eventually,  in  Regime  II,  the  free  electron  concentration  becomes 
such  that  inelastic  electron-atom  collisions  dominate  the  reaction  and  the  degree 
of  ionization  approaches  a  maximum  corresponding  to  Saha  equilibrium  (Ref ..-58). 

The  onset  of  Regime  II-I,  at  a  point  just  prior  to  Saha  equilibrium  is  characterized 
by  a  strong  emission  of  radiation  (Refs.  59  j  60)  which  results  in  plasma  cooling 
and  hence  in  electron-ion  recombination.  This  effect  is  illustrated  in  Fig.  6 
by  the  slow  rise  in  the  fringe  shifts  and  the  corresponding  fall  in  the  plasma 
temperature  and  degree  of  ionization. 

4.2  Experimental  Measurements 

4.2.1  Flow  Disturbances 

An  initial  series  of  shock  experiments  performed  in  pure  argon  were 
consistently  plagued  by  a  complex  system  of  flow  disturbances  of  wavelength  and 
amplitude  sufficient  to  make  accurate  interpretation  of  the  interfe'rogr.ams'n 
extremely  difficult.  Although  present  on  all  the  interferograms  of  this  series, 
idential  wave  geometry  was  not  reproducible  from  run  to  run  for  otherwise  identical 
boundary  conditions.  Typical  interferograms  shown  in  Fig.  7  for  the  ruby  laser 
and  second  harmonic  wavelengths  clearly  do  not  exhibit  the  stable  behaviour 
anticipated  in  Fig.  6(i). 

Inspection  of  Fig.  7 3  especially  the  interferpgram  for  6943  $ 
indicates  the  existence  two  opposing  symmetrically  displaced  transverse  waves  with 
a  region  of  confluence  in  the  centre.  In  addition  the  visible  non-planarity  of 
the  shock  front  appears  to  be  directly  related  to  this  same  systefii  suggesting 
that  the  transverse  waves  are  weak  forward  moving  oblique  shocks.  It  is  interest¬ 
ing  that  no  references  to  any  such  flow  disturbances  were  made  either  by.Alpher 
and  White  (Ref.  2)  or  by  Bershader  et'  al  (Refs,  58,  59*  66)  both  groups  of  which 
performed  space  resolved  optical  interferometry  on  ionized  argon  shocks  (although 
the  interferogram  in  Fig.  1  of  Ref.  69  does  show  the  presence  of  a  relatively  weak 
system  of  waves  propagating  in  Regime  III) . 

It  is  therefore  significant  that  the  only-  difference  between  thev 
conditions  of  the  present  experiment  and  those  of  the  aforementioned  workers 
are  the  shock  tube  cross  sectional  dimensions.  Alpher  and  White  and  Bershader 
et  al.  used  tubes  of  8.25-cm  x  8.25-cm  and  5-cm.x  5-cm  cross  section  vespectively 
whereas  the  dimensions  of  the  present  facility  are  approximately  l8-cm  x  10-cm. 
This  suggests  the  possibility  of  a  resonance  effect  in  which  the  tube  cross 
sectional  dimensions  and  the  ionization  relaxation  zone  length  are  such  as  to 
form  a  cavity  capable  of  supporting  a  system  of  transverse  waves. 

Figure  8  shows  a  similar  pair  of  interferograms  for  a  shock  propa¬ 
gating  into  nitrogen  under  similar  Mach  number  and  initial  pressure  conditions. 

The  exponential  type  dissociation  relaxation  zone  is  clearly  visible  'for  which 
the  equilibrium  degree  of  dissociation  is  about  24$.  It  is  interesting  to  see 


that  the  flow  is  completely  stable  and  devoid  of  disturbances  with  the  shock  front 
exhibiting  a  high  degree  of  planarity  (as  can  be  seen  by  comparison  with  the  ad¬ 
jacent  reference  wire).  This  suggests  that  the  transverse  wave  phenomenon  already 
described  is  specific  to  ionized  shocks  for  which  the  ratio  of  the  shock  tube:  width 
,dr  -height  to  the  relaxation  zone  length  exceeds  some  critical  value.  As  these  waves 
are  not  observed  behind  dissociating  nitrogen  shocks,  it  seems  reasonable  to' .assume 
that  the  argon  plasma  is  providing  a  means  for  wave  amplification.  •  ' - 

A-.  ,/  Normally,  in  the  absence  of  applied  electric  or  magnetic  fields.,  plasma 

'Instabilities  on  a  macroscopic  scale  occur  due  to  the  existence  of  a  temperature  ’  ■  ■ 
.difference  between  the  free  electron  and  heavy  particle  gases  where  the  latter 
4s  -'composed  of  ions  and  neutral  atoms,,.  The  amplification  of  initially  acoustic, 
.dlstpibances  in  a  collision  dominated  plasma  has  been  treated  by  Morse  arid  ■' 

,  .Sigard.'fRef s .  6l,  62) „  '  u  -A  - 

However,  these  authors  chose  to  examine  the  more  usual  case  in  which', 
fihe. electron  temperature  Te  is  initially  much  larger  than  the  heavy  particle-  - 
'4'gttperature  Ta.  Such  conditions  are  typical  of  electrically  generated  plasmas 
.;ifu Which  the  free  electrons  are  initially  excited  to  a  high  temperature  by  the 
application  of  an  electric  field  which  is  then  terminated  allowing  the  two 
temperature  plasma  to  relax  to  equilibrium,  the  kinetic  energy  of  the  free 
electrons  then  being  transmitted  to  the  heavy  particle  gas  through  either  a  > 
random  or  a  coherent  collision  process, 

A--.,  ;  , 

The  model  of  Morse  and  Ingard -proposes  that  in  the  presence  of  'a  •  -  : ‘ 
low-,  frequency  acoustic  source  where  initially  Te/T?  »  1,  the  electron  energy  - 
Can  become  coupled  to  the  acoustic  wave  such  that  a  positive  feedback  occurs 
allowing  the  electronic  energy  to  be  transformed  into  acoustic  energy.  The 
acoustic  waves  then  undergoes  spontaneous  amplification  as  a  result  of  this 
■Coherent  energy  transfer  where  the  plasma  components  all  move  in  phase,  in 
addition,  the  variation  in  the  state  properties  of  the  plasma  components  arq^  •  - 
all  perturbed  to  the  same  degree.  In  "the  absence  of  an  initial  acoustic  per¬ 
turbation  of  the  plasma  the  thermal  energy  of  the  free  electrons  would  be  transf 
ferred  to  the  heavy  particle  gas  by  random  elastic  collisions. 

It  is  now  relevant  to  discuss  this  wave  amplification  process  as  it  ’might 
apply  to  a  thermally  generated  plasma  of  the  type  encountered  in  the  present 
experiments.  The  temperatures  for  the  free  electron  and  the  heavy  particle  gases 
existing  in  Kegimes  I  and  II  behind  the  incident  shock  front,  are  the  converse 
of  those  discussed  by  Ingard  and  Morse  such  that  initially  Te/T&  «  I. 

Figure  9  shews  a  typical  two  temperature  relaxation  profile  for  a  M  =  18 
shock  into  argon  at  p,  =3  mmHg  as  calculated  by  Oettinger  and  BersSader  (Ref.. 

6o).  1 

Without  any  theoretical  justification  it  is  proposed  that,  in  the  , 
case  of  the  thermally  generated  plasma,  the  reverse  process  is  possible  where 
the  coherent  transfer  of  energy  is  now  directed  from  the  heavy  particle  gas  to 
the  free  electrons.  There  is  however  some  experimental  evidence  to  support  this 
thesis  in  that  the  waves  exist  through  the  two  temperature  nature  of  the  plasma 
existing  in  Regimes  I  and  II. 

Calculations  by  Measures  and.  Belozerov  for  strong  shocks  into 
hydrogen  (Ref.  48)  have  shown  that  the  electron  temperature  Te  is  essentially  .  . 
equal  to  the  heavy  particle  temperature  Ta  in  the  post  shock  relaxation  zone <, 

As  shown  in  Fig.  10,  the  temperature  difference  is  virtually  zero  which,  contrasts 
sharply  with  the  large  difference  between  the  two  temperatures  for  the  case  of 
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argon  as  seen  in  Fig.  9- 


This  difference  in  behaviour  between  the  hydrogen  and  argon  plasmas 
arises  due  to  the  much  larger  elastic  energy  transfer  rate  between  the  free  electrons 
and  hydrogen  in  comparison  to  the  corresponding  value  for  argon.  This  is  due  to 
two  effects.  Firstly  (monatomic)  hydrogen  has  l/4o  of  the  mass  of  argon  and  is 
therefore  in  better  mechanical  contact  with  the  electron  gas.  Secondly,  unlike 
that  for  argon,  the  hydrogen- electron  elastic  collision  cross  section  does.-not 
exhibit  the  Ramsauer  effect  (Ref.  63).  This  phenomenon  constitutes  an  anomalous 
minimum  in  the  electron-atom  elastic  collision  cross  sections  for  the  heavier 
noble  gases,  which  in  the  case  of  argon  occurs  at  about  0.4  eV.  However  this 
minimum  is  relatively  broad  so  that  for  a  thermal  distribution  of  velocities  at 
13500°K,  the  transfer  of  kinetic  energy  to  the  electron  gas  from  the  neutral  atom 
gas  through  the  process  of  elastic  collisions  is  relatively  inefficient.  It  is 
therefore  clear  that  the  influence  of  the  Ramsauer  effect  on  the  relaxation  process 
will  be  felt  primarily  in  Regime  II  where  the  temperatures  for  both  the  free 
electrons  and  the  neutral  atoms  are  approaching  the  equilibrium  value  (see  Fig, 
6(ii)). 


It  was  therefore  proposed  to  add  a  small  percentage  of  hydrogen 
to  the  argon  test  gas  with  the  purpose  of  reducing  the  temperature  difference 
existing  between  the  electrons  and  the  .'neutral  atoms  and  ions  sufficiently  to 
create  conditions  unfavourable  for  wave  amolification.  The  hydrogen  additive, 
which  for  the  present  shock  conditions  becomes  completely  dissociated  oh  passing 
through  the  translational  shock,  has  the  effect  of  increasing  the  rate  of . energy 
transfer  from  the  argon  atoms  to  the  free  electrons  through  a  random  elastic 
collision  process  thereby  raising  the  free  electron  temperature  to  a  value  closer 
to  that  for  the  neutral  atom  gas.  Simultaneously  this  should  tend  to  eliminate 
conditions  suitable  for  the  acoustic  wave  energy  transfer  mechanism  which  appears 
to  be  highly  competitive  in  the  case  of  pure  argon. 

A  series  of  experiments  was  therefore  conducted  in  which  varying 
percentages  of  hydrogen  were  added  to  the  argon  test  gas  for  constant  initial 
conditions  of  M  =17-10  and  p  =  2.85  mmHg.  Typical  argon  shock  interferograms 
are  shown  in  Fig.  11  and  Fig.  12  corresponding  to  concentrations  by  pressure  of 
0.033}  0.2  and  0.4$  of  molecular  hydrogen  in  the  pre-shock  test  gas.  The  improved 
planarity  of  the  shock  front  and  stability  of  the  relaxation  zone  with  increasing 
hydrogen  concentration  is  clearly  visible. 

Several  interesting  effects  were  observed.  Firstly  ir  'ing 
the  percentage  of  hydrogen  in  the  test  gas  had  the  desired  effect  of  damping 
out  the  previously  observed  wave  disturbances  until  conplete  stability  of  the 
shock  relaxation  zone  was  achieved  for  a  hydrogen  concentration  of  0.4$  by 
pressure. 


Secondly  as  anticipated  the  relaxation  length  for  Regime  I  becomes 
considerably  reduced  as  the  atom-atom  ionization  cross  section-  for  hydrogen  is 
much  larger  than  that  for  argon.  Specifically  the  initial  slope  of  the  hydrogen- 
hydrogen  excitation  collision  cross  section  curve  for  the  rate  controlling  ground 
state  to  first  excited  state  transition  (Ref.  48)  is  about  two  orders  of  magnitude 
larger  than  the  corresponding  value  for  argon  (Ref.  64).  As  already  postulated, 
the  presence  of  small  percentages  of  hydrogen  in  the  argon  test  gas  will  raise 
the  electron  temperature  both  in  Regimes  I  and  II.  It  can  therefore  be  expected 
that  the  lengths  of  Regime  II  will  also  be  reduced  as  the  ionization  rate-  for 
electron-atom  collisions  is  strongly  dependent  on  the  electron  temperature. 
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varying  approximately  as  Te'  (Ref.  58). 

Figure  13  shows  a  plot  of  the  normalized  total  relaxation  zone 
length  p.X  for  argon  shocks  as  a  function  of  the  concentration  of  molecular 
hydrogen  in  the  initial  test  gas  for  M  =  17  and  p.  =  2.85  mmHg  where  X  is  the 
experimentally  measured  overall  relaxation  zone  length.  The  (visually  fitted) 
curve  tnrough  the  data  points  appears  to  asymptote  satisfactorily  towards  the 
value  for  pure  argon  (10  ppm  inpurity  level)  obtained  by  Wong  and  Bershadep  (Ref. 
58) -»  The  addition  of  relatively  large  concentrations  of  oxygen  or  helium  had  little 
effect  either  on  the  wave  disturbances  or  the  relaxation  zone  length  (see  Fig. 13)* 

A  further  series  of  shock  runs  indicated  that  the  addition  of  0 ,4$ 
of  molecular  hydrogen  to  the  argon  test  gas  resulted  in  the  complete  damping  of 
the  previously  observed  disturbances  for  Mg  =  17  to  24  and  for  corresponding 
p^  =  3.0  to  0.3  mmHg.  Consequently  all  shock  runs  made  for  the  purpose  of  measur¬ 
ing  the  refractivity  of  singly  ionized  argon,  were  made  using  a  test  gas  contain¬ 
ing  0.-4$  of  molecular  hydrogen  by  partial  pressure. 

The  enployment  of  this  expedient  immediately  raises  the  question 
concerning  the  effect  of  the  hydrogen  additive  on  the  experimental  fringe  shift 
data,  on  both  the  experimental  and  calculated  shock  properties  used  in  the  data 
reduction  and  consequently  on  the  final  refractivity  results  for  singly  ionized 
argon.  A  calculation  was  therefore  performed  to  estimate  the.  effect  pf  0.4$  of 
hydrogen  on  the  frozen  argon  shock  properties  as  exist  immediately  behind  the 
translational  shock  front  prior  to  the  onset  of  ionization.  It  was  assumed  that 
the  molecular  hydrogen  becomes  completely  dissociated  on  passigg  through  the  inci¬ 
dent  shock.  For  shock  Mach  numbers  from  17  to  2k  and  for  the  corresponding  frozen 
shock  tenperatures  ranging  from  25,000  to  50,000  K,  the  assumption  of  rapid,  and. 
complete  dissociation  of  the  small  hydrogen  concentration  would  seem- plausible. 

The  results  of  this  calculation  are  presented  in  Fig.  14  as  percentage  deviations 
from  the  pure  argon  case.  Changes  in  the  post  shock  pressure  are  negligible 
whereas  the  values  for  the  density  and  tenperature  are  approximately  +  0 .4$  and 
-0,7%  respectively  for  M  =  17.  These  deviations  become  smaller  with  increasing 
Mach  number  primarily  because  the  dissociation  energy  for  hydrogen  becomes  a,., 
progressively  smaller  fraction  of  the  total  energy..  A  similar  calculation  for 
Saha  equilibrium  was  not  performed  as  it  was  considered  that  the.  results  would 
reflect  the  figures  already  obtained  for  the  frozen  shock  state,  particularly 
as  the  ionization  potential  for  hydrogen  (13.6  eV)  is  of  a •  similar  magnitude  to 
that  for  argon  (15.8  eV). 

Consequently  as  the  effects  of  the  0.4$  of  hydrogen  on  the  frozen 
shock  properties  for  argon  are  small,  it  is  suggested  that  the  errors  incurred 
by  ignoring  the  presence  of  the  hydrogen  on  the  shock  properties  and  the  equili¬ 
brium  fringe  shift  measurements  would  be  considerably  smaller  than  those  from 
other  sources  of  error.  The  analysis  of  the  refractivity  data  was  therefore 
performed  on  the  assunption  that  the  -ringe  shift  measurements  were  made  for  pure 
argon  using  equilibrium  shock  Hugoniot  data  calculated  for  the  case  of  pure  argon. 

So  far  no  comment  has  been  made  concerning  either  the  source  of 
the  initial  perturbation  for  the  flow  disturbances  or  the  wave  structure  of  the 
unsteady  ionized  shock  flow.  There  exist  several  types  of  initial  perturbation 
sources  capable  of  generating  acoustic  waves.  These  can  conveniently  be  Bub*r 
divided  into  three  categories  as  each  type  tends  to  act  on  the  flow  in  a  different 
way. 


The  first  group  consists  of  stationary 'sourc.es  located  near  the  test 
section  such  as  small  wall  imperfections,  window  bevels,  pressure  transducer- 
plugs,  etc.,  which  are  capable  of  producing' Mach  waves  in.  the-  supersonic. -post  - 
shock  flow.  Although  much  effort  was  expended  in  reducing  the*  number  and  magni¬ 
tude  of  these  wall  imperfections,  this  had  ho  apparent  effect  in  improving-  the 
stability  of  the  shock  flow.  In  addition  it  is  difficult  to  visualize  how  the 
observed  disturbances  as  typified  by  the  interferograms  in  Figs’.  7  and  11  might 
develop  in  the  region  of  the  test  section  during  the  period  of  shock  passage, 
particularly  in  view  of  the  stable  nature  of  the  shock  flow  in  the  case  of  the 
dissociating  nitrogen  shock  shown  in  Fig.  8.,  The1  second1  group  of  perturbations 
can  be  considered  as  acoustic  or  weak  shock  disturbances  produced. -either  at -or 
upstream  of  the  contact  surface.  The  contact  surface  is  generally  a  diffuse- 
turbulent  region  which  might  constitute  a  weak  acoustic  source.,  -In  addition 
there  exists  the  possibility  of  a  Rayleigh-T’aylor  contact  surface  -instability 
(Ref.  65)  which  may  be  specific  to  combustion  driving  into1  the  heavier  argon 
test  gas  but  not  into  the  lighter  nitrogen- test.  gas.  Combustion  instabilities 
in  the  driver  gas  particularly  in  the  case  of  constant  pressure  combustion  where 
combustion  continues  after  diaphragm  rupture,  has  been  Observed,  to-  generate 
unsteady  flow  in  the  driven  shocked  gas  (Ref;  66);;  The  third'  and  final-  group 
of  flow  perturbations  are  those  classified  as  remote  stationary  sources  located 
at  distances  greater  than  60  (equi valent)  tube  diameters-  upstream  of  the  test 
section.  Such  features  are  the  two  2-in.  dia.  vacuum  pumping  ports-,,  the- round 
to  rectangular  transition  section  and  -the  ruptured  diaphragm.  Disturbance  from 
sources  of  this  type  generally  propagate -downstream  via  a  series  of  multiple 
transverse  Mach  wave  reflections  (-Ref.  67 )  which-  in'  the  -absence  of  any  mechanism 
for  amplification,  decay  -as  x~’3/2,  where  x  is  the'  'distance-  from  the  source  to-  - 
the  point  of  observation  (Ref.  68)'-.  •  " 

In  the  case  of  the  ionized  shock  flow,'  it  is -possible "for  these  per¬ 
turbation  sources  to  act  in  two  ways.  Firstly,  the  flow  can  !pick-upr  th§ 
disturbances  as-  it  passes  the  perturbation,  source  which  is  then-  sustained  at  a  -  '• 
given  level  by  wave  amplification  as  the  shock  flow  progresses  down  the  tube. 
Alternatively,  the  disturbance  can  be  transmitted  through  the  flow  in  a  manner 
similar  to  the  transverse  Mach  waves  referred  to  above.' 

Allithree  types  of  perturbation  source,  with  one'  possible  exception 
should  in  principle  generate  a  system  of  disturbances  in.  the- -form  of  either 
transverse  or  normal  waves  which  should  "be  detectable  in  any- equilibrium  shock- 
flow.  It  is  therefore  significant  that  wave  structure  was  hot  observed  in  the 
nitrogen  shock  flow  interferograms  shown  in  Fig.-  8.  However  it  must  be  recognized 
that  a  schlieren  picture  of  the  nitrogen  flow  -might  reveal  the  existence  of "wave 
structure  not  visible  on  an  ihterferbgram*.  The- exception  referred  to  is,  the- ease' 
in  which  the  initial  disturbance  is  generated  at  a.  point -remote  from-  the  test 
section.  This  system  is  then  amplified  by  the  ionized  shock  flow  ’as  it  is  swept 
downstream  to  the  test  section  such  that  a  stable  system,  of  transverse  waves  is 
maintained  where  all  knowledge  of  Original  source  is  lost;  This  latter  mechanism 
for  wave  initiation  would  seem  to  be  the  most  plausible  in  ’view  Of  the -existing 
experimental  evidence.  But  clearly  a  carefully  controlled  experimental  investi¬ 
gation  of  this  phenomenon  is  required  in  order  to  locate  the ’-origin  of  -this  -un¬ 
steady  shock  flow;  The  interferograms  presented  in  this-  work  do  however  contain 
information  regarding  the  mode  of  propagation  of  the  disturbances  in  the  shock 
front  and  various  ionization  regions.'  ' 

*  Figure  8(a)'  shows  a  similar  interfer.ogram.- for- a  :M  -= -11.6  nhock  into  oxygen 
at  p  =  bO.O  mmHg  for  which  the  equilibrium- degree  of  dissociation,  is  .about  18^. 
Due  zo  the  relatively  high  initial  pressure,  the,  disturbances  referred  to  are 
now  visible  as  a  system  of  transverse  waves.-. 
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If  it  assumed,  that  the  initial  disturbance  is  either  an  acoustic  or 
a  weak  compression  wave,  propagating  towards  the  shock  front,  then  according  to 
the  arguments  presented,  amplification  should  not  occur  until  the  disturbance 
enters  Regime  II  (see  Fig.  6(iii)-).  Although  Regime  III  is  one  in  .which  piasma 
cooling  occurs  due  to  radiative  losses  leading  to  recombination,  there  appears 
to  be  no  reason  to  expect  the  election  temperature  to  differ  from  that  for  the 
heavy  particles,  suggesting  therefore  that  wave  amplification  does  not  occur  in 
this  flow  regime. 

It  would  now  seem  reasonable  to  suggest  that,  wave  amplification  varies 
at  some  function  of  the  product  of  the  temperature  difference  (Ta-Te)  and  the 
electron  density  We  so  that  wave  amplification  will  attain  a  maximum  at  some  point 
in  Regime  II  where  We  is  rising  concurrently  with  the  fall  in  (Ta-Te)  as  equilibrium 
is  approached.  Inspection  of  the  interferograms  in  Fig.  7  and  ll(i)  tends  to 
support  this  thesis  in  so  much  as  the  disturbances  appear  to  achieye  their  maximum 
amplitude  in  this  region  particularly  at  the  poi-njbs  of  wave  confluence  located  be¬ 
hind  the  shock  front  wave  heads..  The  fact  that  incident  sho.ck  front  has.  become 
locally  intensified  by  wave  interaction  indicates  that,  the  disturbances  .have  be¬ 
come  weak  shock  waves  (Ref,.  ,69),  presumably  due  to  wave  amplification. 

Under  these  conditions,  the  interaction , also  produces,  a  .rearward 
moving  rarefaction  wave  (Ref.  69),  which  then  interacts  with  and, intensifies  both 
the  weak  incoming  shocks  and  the  ionization  density  jump  whilst  the  rarefaction 
wave  itself  becomes  attenuated  (Ref.  69).,  It  is  possible  that  these  quasi-normal 
waves  visible  in  Regime  III  of  the  interferograms  in  iFig..7  and  11,]  are  in  fact 
tne  rearward  moving  rarefaction  waves  referred  to  above.  Why  these  waves  should 
decay  so  rapidly  as  they  progress  upstream  into  Regime  ±11  is  not  .clear  unless  wave 
damping  is  accomplished  through  a  coupling  with  the  radiative  cooling  process,. 
However,  it  must  be  pointed  out .that  it  is  .not  proven  that  these  waves  are  in  fact 
rearward  moving;  multiple  frame  interferometry  of  the  .shock  flow  would  Vesolve  this 
question.  (W.B.  A  similar -group- of  waves  is  visible  .in  Regime  III  of  the  ionized 
argon  shock  interferogram.  shown  in.. 'Fig..  1  of  Ref .  66).  _ 

There  is  however  one  aspect  of  the  observed  flow  disturbances ,  so  far. , 
not  discussed,  which  lends  some  creadence  to. the  foregoing  qualitative  analysis , 

It  was  noticed  that  there  was  a  systematic  tendency  for  the  number,  of  shock  induced 
heads  on^the  incident  shock  front,  to  increase  whilst  simultaneously  becoming  smaller 
in  amplitude  with  reduction  in  the  ionization  relaxation  zone,  length,  due  to. the 
effect  of  the  hydrogen  additive  on  the  ionization  process.  The  shocks  shown  in 
Fig.  7,  ll(i)  and  ll(ii)  all  for  similar  boundary  conditions  except  _for  the  hydro¬ 
gen  concentrations  of  0,  p. 033  and  0.4$.  respectively,,  have  1,  3  and *4  heads  res¬ 
pectively.  It  therefore  appears  that  the  region. bounded  by  the  shock  front  and 
the  ionization-induced  density  rise  is  tending  to  behave  as  ..a  resonant  cavity  in 
which  the  initial  wave  perturbation  becomes  intensified  both,  by  the  acoustic  wave 
amplification  process  and  -by  coalescence  with  additional  waves  entering  this  zone. 
This  suggests  that  the  ratio  of  the  shock  tube  width. to  the  ionization,  relaxation 
zone  length  in  conjunction  with  the  Mach  wave  angle  is  a  .governing  factor  in  this 
process,  particularly  as  these  waves  were- not  observed  in  the  argon  shock  flow 
interferograms  of  Bershader  et  a-1  (Refs.  58,  ,!?9>  6Q\  and  of  Alpher  and  White 
(Ref.  2).  .  ...  "  .  ’ 

T.j  is  interesting  that  this  behaviour,  particularly  as  exemplified 
by  the  shock  flow  in  Fig.  ll(i).  is  similar  to.  that,  observed  by  .White.. (Ref,.,  ,7<D). 
for  detonations  ip  stoichiometric  hydrogen- oxygen  mixtures-  diluted  in  xenon,  A 
resonant  transverse  wave' structure  was  also-  obtained  for  which  the  number  of  wave 
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heads  was  considered  to  be  related  to  the  reaction  zone  thickness. 

It  is  realized  that  the  foregoing  qualitative  discussion  pertaining 
to  the  origin  and  behaviour  of  these  flow  disturbances  poses  as  many  questions  as 
it  answers  and  must  therefore  only  be  considered  as  tentative.  Clearly  then  an 
experimental  and  a  theoretical  investigation  of  this  phenomenon  must  be  conducted 
with  particular  emphasis  being  put  on  the  following  points: 

a)  Conduct  an  experimental  investigation  to  locate  the  initial  perturbation 
source. 

b)  Perform  a  theoretical  study  of  the  acoustic  wave  amplification. process  for 
collision  dominated  plasmas  where  (Ta/Te)»  1  and  N  =  10^7  cm” 3. 

c)  Conduct  an  experimental  and  theoretical  examination  of  the  action  of  the 
hydrogen  additive. 

(i)  on  the  nonequilibrium  temperature  difference, 

(ii)  on  the  kinetics  of  regimes  I,  II  and  III  and 

(iii)  on  the  acoustic  wave  anplification  process. . 

and 

d)  Analyze  both  theoretically  and  experimentally  the  various  wave  interaction 
pr-  cesses  which  have  been  observed  to  exist  in  the  present  ionized  argon 
shock  flows. 

4.2.2  Interferometric  Measurements 

A  total  of  ,23  shock  funs  were  performed  in  argon  containing  0 .4$ 
of  molecular  hydrogen  by  partial  pressure  with  flow  and  no-flow  interferograms 
obtained  at  the  wavelengths  of  the  ruby  laser  (6943a)  and  its  second  harmonic 
(3471a).  The  shock  boundary  conditions  ranged  from  M  =  17  to  24  arid  p.  =  3 
to  0.3  mmHg  with  corresponding  degrees  of  ionization  from  17$  to  48$, electron 
densities  £>f  about  101'  cm”3  arid  pressures  of  about  1  atmosphere.  The  equilibrium 
plasma  temperature  remained  virtually  constant  at  about  i.3500°K  for  ail  the 
shock  runs  due  to  the  compensatory  effect  of  the  rise  in  the  degree  of  ionization 
with  increasing  shock  Mach  number. 

The  fringe  shifts  between  preshock  stat’e  1  and  Saha  equilibrium 

state  2  were  obtained  according  to  the  procedure  outlined  in  Appendix  A  using 
Eq.  A. 4.  The  location  of  the  point  corresponding  to  Saha  equilibrium  along  any 
given  fringe  was  taken  to  be  that  at  which  the  plasma  refractive  index  for 
state  2  attains  a  minimum.  With  reference  to  the  interferograms  in  Figs..  12 
and  15,  Saha  equilibrium  was  taken  to  exist  when  the  negative  fringe  S'V,  reaches 

it®  maximum  numerical  value  prior  to  the  onset  of  radiation  cooling.  Values  of 
Si'p  were  obtained  for  four  vertically  distributed  locations  across  the  shock  .wave 
interferograms  for  a  given  wavelength  and  shock  run.  Hence  as  the  valiie  of  *12 
given  by  Eq.  A. 4  is  itself  the  average  of  two  measurements,  then  the  final  value 
is  the  arithmetic  mean  of  eight  individual  measurements.. 

As  dis-unsed  in  Section  31.3 ,  it  was  considered  possible  to  locate 
the  fringe  centre  lines  ;o  an  accuracy  of  +  0.2  mm  which  corresponds  to  a  maximum 
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uncertainty  in  the  fringe  shift  measurements  of  +  0.04  fringes  for  X  =  3471S 
for  which  the  fringe  spacing  is  a  minimum. 

An  additional  seven  runs  were  performed  using  the  neodymium  glass 
laser  (10,600a)  and  its  second  harmonic  generator  (5>300$)  thereby  utilizing  the 
strong  wavelength  dependence  of  the  free  electron  refractivity  (see  Eq,  2.1.5)  in 
order  to  obtain  accurate  experimental  measurements  of  the  electron  density. 


Typical  shock  flow  interferograms  for  the  ruby  and  neodymium  lasers 
for  this  final  series  are  shown  in  Figs.  12  and  15  respectively.  The  large  equili¬ 
brium  fringe  shift  (S^?  =  7)  for  the  neodymium  laser  wavelength  (10,600$)  and  the 
marked  effect  due  to  tne  radiation  cooling  in  Regime  III  demonstrates  the  large 
negative  contribution  of  the  free  electrons  to  the  plasma  refractive  index.  The 
corresponding  (positive)  fringe  shift  across  the  frozen  shock  (S^  =  0.2)  is  due 
solely  to  the  compression  of  the  neutral  atoms. 

The  analysis  and  discussion  of  the  electron  density  and  ion  refrac¬ 
tivity  measurements  are  given  in  Sections  4.3  and  4.4  respectively. 

4,3  The  Experimental  Determination  of  the  Equilibrium  Electron  Density 

Accurate  values  for  the  neutral  atoms,  first  ion  and  free  electron 
densities  are  clearly  necessary  if  the  refractivity  for  (singly)  ionized  argon  is 
to  be  measured  by  the  interferometric  method  described  in  Section  2.2.  The  experi¬ 
mentally  determined  free  electron  density  is  obtained  by  the  two-wavelength  inter¬ 
ferometric  method  described  in  Section  2.3  which  is  independent  of  the  existence 
of  either  complete  or  local  thermodynamic  equilibrium  fbr  the  plasma.  Hence  this 
measured  value  can  be  used  as  an  indicator  to  show  whether  the  concentrations  of 
the  plasma  constituents  have  achieved  the  values  given  by  the  Saha  equation  for 
the  given  boundary  conditions.  As  electron- ion  recombination  due  to  radiative 
cooling  is  known  to  commence  at  some  point  in  Regime  II  (see  Fig.  6)  prior  to 
equilibrium  (Ref.  60) ,  it  is  possible  for  the  measured  equilibrium  electron  density 
to  lie  below  the  Saha,  equilibrium  value  at  all  times  iii  which  case  it  would  not  be 
possible  to  determine  the  ion  refractivity  from  Eq.  2.2,4. 


For  the  laser  fundamental  wavelength.  Ay  and  its  corresponding  second 
harmonic  Ax  (=  Ay/2),  NE  the  experimental,  electron  density  is  obtained  from 
Eq.  2.3.2  as  c" 


-  N  A 
=-2_X 

e  A_ 


(p  y-l)  L 
e  ' o 


(M-l) 


where  it  is  assumed  that  the  ratio  of  the  ion  refractivity  to  the  atom  refractivity 
is  independent  of  wavelength.  The  experimental  fringe  shift  terms,  sV  and  S^  for 

the  second  harmonic  and  fundamental  respectively  were  obtained  as  described  in. 
Section  4.2.2  and  Appendix  A.  L  is  the  test  section  geometrical  path  length  and 
is  equal  to  4.010  inches. 
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The  experimental  electron  densities  were  compared  to  two  theoretical 
values,  both  of  which  were  calculated  for  Saha  equilibrium  assuming  radiation  losses 
to  be  negligible  and  using  the  shock  Hugoniot  equations  for  the  given  (experimental) 
values  of  M  ,  p,  and  T  .  The  first  calculated  value,  R,  G,  was  obtained  using  , 

electronic  excitation  partition  functions  which  included  only  the  ground  state 
terms.  The  second  and  more  realistic  calculated  value,  N  ^  was  obtained  using 

G 

partition  functions  which  included  terms  for  the  excited  electronic  states  in  addi¬ 
tion  to  that  for  the  ground  state.  The  effects  of  the  plasma  microfield  on  the 
neutral  atom  partition  function  and  the  ionization  potential  were  estimated  accord¬ 
ing  to  the  method  of  Margenau  and  Lewis  (Refs,  jl,  72)  in  which  ail  those  electronic 
orbits  were  suppressed  whose  classical  major  semi-axes  exceeded  the  Debye  shielding 
distance.  Data  for  the  energy  levels  of  neutral  and  singly  ionized  argon  up  to  a 
principal  quantum  number  of  32  were  taken  from  the  tabulations  of  Moore  (Ref.  25) 
together  with  predicted  values  taken  from  Ref.  73  • 


Before  comparing  R  to  these  calculated  values,  it  is  pertinent  to 
discuss  how  and  why  these  two  calculations  differ  from  each  .othfer...  N-  * .  and  N  G  are 

compared  on  a  percentage  basis  in  Fig.  16  in  which  (R.  -R  ;/N  r  is  plotted- as  a 
function  of  M  for  p  =  10,  1  and  0.1  mmHg.  Several,  f oirits  ofein'ter'est  apparent 
from  this  graph  are  aiscussed  below. 


G  F  G 

(i)  Nu  is  always  greater  than  Ng  because  in  the  calculation  of  Rg  ,  the 

omission  of  the  excited  electronic  state, ..terms  from  the-  partition 
functions  for  the  atom  and  ion  results  in  an  increase  in  energy  avail¬ 
able  for  investment  in  ionization.  In  addition  the  ratio,  of  the  ion 
partition  function  to  that  for  the  atom  as.  psed  in  the  Saha  equation 
increases  due  to  the  relatively  large  excitation  energies  for  the 
excited  states  of  singly  ionized  argon,  resulting  in  ah  increase  in 

R  G  in  relation  to  R  F. 
e  e’ 

G  F 

(ii)  (Wg  -  Ne  )  increases  with  p^  for  a  given  because  as  the  degree  of 

ionization  decreases  with  increasing  p^,  more  energy  becomes  available 
for  populating  the  excited  electronic  states  in  the  case  of  as  this 

is  a  non-pressure  sensitive  process  (if  to  a  first  approximation  one 
ignores  the  plasma  microfield  effect).  Since  this  effect  does  not 
exist  in  the  case  of  -  ,  then  (RgG-  N^-  )  increases  with  p.  for  a 

given  value  of  Mg.  e 

G  F 

(iii)  (R  -  R  )  increases  with  M  for  a  given  p. .  This  can  be  understood 

G  G  S  X 

from  (i)  above,  in  which  thG  excited  state  population  in  the  case  .of 
R  ^  increases  with  increasing  M  .  •  „ , „  -  -  . 

G  F  F 

(iv)  It  is  observed  from  Fig.  l6  that  all  of  the  curves  of  ('Nv7-  ’N-  )/Ng 
plotted  for  constant  p^  attain  a  maximum  at  a  given  value -of  Mgi 

This  presumably  occurs  because  as  the  neutral  atom  concentration 
falls  with  increasing  Mg  for  constant  p^,  then-  the  total  energy 
invested  in  populating  the  excited  states,  of  neutral  argon'  (in  the 
case  of  N  F)  also  falls  thereby  reducing  (R  G  -  R1')-.  Due  to  the 

large  excitation  energies  for  singly  ionized  argon  in  relation  to 
those  for  neutral  argon,  the  excited  states  for  this  species  do  not 
become  commensurately  populated  as 'the  degree  of  ionization  increases. 
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G  F  F 

Values  of  (N  -  N  )/N  for  the  thirty  experimental  shock  runs 
6  ®  ®  G  F  F 

are  shown  in  Fig.  17  as  a  function  as  Mg.  The  average  value  of  (Ng  -  N  )/N£ 

is  approximately  +  3.5$  with  the  curve  exhibiting  a  slow  decrease  with  increase 
in  M  due  to  a  corresponding  drop  in  p  (for  the  reasons  given  in  (ii)  above). 

S  X 


The  principal  ourp os e  of  Fig.  17  is  to  compare  the  experimentally 
measured  electron  density  nJ* to  the  calculated  value  N^,  On  average  N  ®  =  1.02. 

N  F)/NeF  =  i.l.9$  where  the  latter  difference 

expression  exhibits  a  slow  fall  with  increasing  M  for  the  reasons  given  in  (ii) 
above. 


F  ^  F 

N  *  +  3$  or  more  specifically  (N 


TjJ 

Discussion  of  this  discrepancy  between  and  Ng  is  treated  in 

two  parts  with  consideration  being  given  first  to  possible  sources  of  systematic 
error  in  and  then  secondly  to  those  in  Ng 

F 

Sources  of  Error  in  Calculating  N_ 


E 


a)  An  obvious  source  of  error  in  N  *  exists  in  model  used  to  terminate  the 
electronic  excitation  partition  function.  For  an  argon  plasma  at  33  ,500  K 
and  1  atmosphere  pressure,  the  model  of  Margenau  and  Lewis  (Refs.  71,  72) 
terminates  the  neutral  atom  partition  function  at  a  principal  quantum 
number  q  =  21.  Alternatively  if  one  reduces  the  ionization  potential 
according  to  the  Debye-Huckel  method  (Ref.  7*0  and  then  terminates  the 
partition  function  at  this  new  ionization  limit,  a  principal  quantum  ^ 
number  q  =  9  is  obtained.  Clearly  this  will  produce  a  reduction  in  N 

by  increasing  the  energy  available  for  ionization  thereby  bringing  N  F 
into  closer  agreement  with  nJO 

F 

b)  As  N  is  related  to  a  particular  experiment  through  the  measured  values 

e  p 

of  M  ,  p,  and  Tn ,  then  errors  in  these  values  generate  errors  in  N 
si  1  E  ® 

in  relation  to  W  which  is  a  product  of  the  fringe  shift  measuiements 
As  indicated  in  eSection  3-1,  the  maximum  errors  in  measuring^Ms  and  p^ 
are  +  1$  and  +  2-l/2$  respectively  which  produce  errors  in  Ne  of  about 
2-l/2$  for  conditions  typical  of  these  experiments.  However  these  error' 
limits  on  N  F  are  due  to  random  errors  and  cannot  therefore  be  considered 
to  be  uhe  source  of  the  systematic  1.9$  difference  between  and  N  ^ 

although  they  are  likely  the  major  source  of  the  +  3$  scatter  in  the 


(V5 


data  in  Fig  .-17. 


Sources  of  Error  Incurred  in  the  Measurement  of  N 


E 


TTjr- 


E 


a)  Ne~  was  obtained  ^rom  the  experimental  fringe  shift  data  on  the  assumption 
that  the  ratio  (pj  -l)  -l)Q  Is  independent  of  wavelength.  As  the 

values  of  N  ®  obtained  using  both  the  ruby  and  neodymium  glass  lasers  are 
consistent,  it  is  doubtful  whether  this  assumption  contributes  to  this 
discrepancy.  Additional  confirmation  of  this  conclusion  (as  wiil  be  shown 
in^the  following  section)  is  provided  by  the  fact  that  the  ratip 
(Pj  -l)o/(p^  -l)Q  is  effectively  constant  in  the  optical  and  near  infrared 
regions  of  the  spectrum. 
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b) 


It  was  shown  in  Section  2.1,  that  the  excited  states  make  a  -6.6$  contribution 
to  the  total  neutral  argon  refractivity  at  6943X  for  conditions  of  13j500°K 
and  1  atmosphere  pressure.  If  this  excited  state  contribution  is  introduced 
into  the  data  analysis  for  the.. ruby  laser  experiments,  it  produces  an  approxi¬ 
mately  1 $  reduction  in  N  E  thereby  halving  the  observed  discrepancy  between 
E  F  ® 

N  and  .  However  the  same  excited  state  refractivity  calculation  for 
neutral  argon  indicated  a  negligible  contribution  at  10,60oX.  Consequently 
as  the  average  valuesoof (NgE  -  NeE)/NgF  for  wavelengths  of  6943$  apd  10,600a 

are  +  1.89$  and  +  2.22 $  respectively  where  this  small  difference  is  in, the 
direction  in-consistent  with  the  calculated  excited  state  refractivity  con¬ 
tributions,  it  is  concluded  that  the  excited^states  are  not  the  principal 

cause  of  the  observed  discrepancy  between  N  **  and  N  * . 

e  e 

E 

The  prime  purpose  for  monitoring  N  was  to  observe  whether  there  was  any 
significant  lowering  of  the  electrSn  density  from. the  value  predicted  for 
Saha  equilibrium  due  to  radiative  recombination.  The  intensity  of  the 
emitted  radiation  is  known  to  follow  the  growth  of  the  free  electron  con¬ 
centration  in  Regime  II,  attaining  a  maximum  at  the  point  corresponding  to 
Saha  equilibrium  (Ref.  60).  From  thereon  the  loss  of  energy  through  radia¬ 
tive  cooling  induces  a  fall  in  the  electron  density;  this  is  characterized 
by  the  relatively  slow  rise  of  the  fringes  in  Regime  III  of  the  interfero- 
grams  in  Fig.  12  and  15. 


E 

It  was  therefore  interesting  to  find  that  the  experimental  values,  N 

b 1  G 

were  consistently  larger  than  N  although  always  lower  than  the  upper  . 
bound  represented  by  N  i  NgE  e  for  the  thirty  shock  runs  of  this 
experiment  are  shown  in  Fig.  l8  as  a  function  of  'Mg  together  'With  N  F 

for  p.  =  10,  1  and  0.1  mmHg.  In  all  cases  N  E  lies  below  the  value  of 
1+17  ® 

7  x  10  cm  J  suggested  by  McWhirter  (Ref.  75)  to  be  necessary  for  the 

attainment  of  local  tnermodynaraic  equilibrium  (LTE)  in  an  argon  plasma  - 
at  an  electron  temperature  of  13,500°K.  Even  if  LTE  does  not  exist  due 
to  radiative  transitions  creating  a  relative  overpopulation  of  the  ground 
state  with  respect  to  a  Boltzmann  distribution,  it  is  clear  that  the 
total  neutral  argon  population  must  be  in  good-  agreement  with  the  Saha 
value  due  to  the  close  correspondence  between 'NgE  and  Ng  .  Hence  any 
deviations  from  LTE  will  increase  the  population  of  the  ground  state  in 
relation  to  the  overall  excited  state  population  thereby  reducing  the 
influence  of  the  large  refractivity  terms  for  these  states  on  the. re¬ 
fractivity  measurements  for  singly  ionized  argon.  However  inspection 
of  the  argon  shock  flow  interferograms  in  Fig.  12  and  15  shows  that  the 
collision  induced  electron-atom  ionization  relaxation  process  of  Regime 
II  is  very  much  more  rapid  than'  the  radiative  recombination  process  in 
Regime  III  due  in  part  to  the  presence  of  the  hydrogen  additive.  Jt'  would 
therefore  seem  plausible  to  suggest  that  under  these  conditions,  the  plasma 
achieves  a  Boltzmann  distribution  of  internal  states  at  the  point 
corresponding  to  Saha  equilibrium  which  is  then  maintained  as  the  plasma  _ 
cools  on  passage  through  Regime  III. 


d)  Griem  (Ref.  76,  see  also  Refs.  13,  l4)  has  suggested- that  electrons  bound 

in  orbits  lying  close  to  the  ionization  continuum,  behave  as  if  free  in 

the  presence  of  electro-magnetic  radiation.  The  effect  of  this  phenomenon 

would  therefore  be  such  as  to  create  an  apparent  increase  in  the  measured 

value  of  N  . 

e 


e)  As  mentioned  in  Section  3»3>  the  fringe  shifts  at  10,60oS  and  6943I  were 
measured  to  an  estimated  accuracy  of  +  0.02  fringes.  For  these  measure¬ 
ments  at  10,600a  and  6943a,  the  minimum  fringe  shifts  were  of  the  order  of 
5  and  2  respectively,  which  produce  maximum  errors  in  N  ®  of  +  l/2$  and 
+  1 %  respectively.  This  is  confirmed  by  the  reduced  scltter  of  the 

(#  -  #  )/M^  data  obtained  using  the  neodymium  glass  laser  in  relation  to 

that  for  the  ruby  laser.  Although  these  error  contribute  to  the  scatter  of 
the  (N  E-  N  E)/n  E  data  as  shown  in  Fig.  17,  they  are  of  a  random  nature 

and  therefore  do  not  contribute  to  the  systematic  difference  between  N 

and  N  F.  e 

e 

E  F  It  is  therefore  suggested  that  the  observed  difference  between 

N  a  and  N  is  due  either  to  the  use  of  an  unrealistic  partition-function  sum  termi¬ 
nation  model  to  calculate  N  r  or  to  the  effect  of  bound  electrons  lying  close  to 
the  continuum  acting  as  if  free.  However  it  is  felt  that  the  experimental  values 
N  E  are  a  better  description  of  reality  in  the  present  experiment,  particularly 
as  these  values  are  obtained  directly  from  the  measured  fringe  shifts.  These  values 
are  therefore  used  in  the  following  section  to  calculate  the  ion  refractivity  term 
(|ij  -l)Q  from  the  experimental  fringe  shift  data. 

4.4  The  Experimental  Determination  of  the  Refractivity  for  Singly  Ionized  Argon 

As  outlined  in  Section  2.2,  the  ion  refractivity  at  STP  is  related 
to  the  experimentally  measured  fringe  shifts  at  Saha  equilibrium  by  Eq.  2.2.4 
which  is 

4  -bo  '  4  -bo  -  (“A  -b0  {(N0/Be,^;p2-P3)/P0)-1  ) 

(M.1) 

where  x,  the  degree  of  single  ionization  has  been  eliminated  through  use  of  the 
equality  x^  =  (We  2/Nq) (p Jp^)  *  T^e  neutral  atom  and  free  electron  refractivities 

at  STP  are  given  by  Eq.  2.1.8  and  2.1.5  respectively.  N  and  p  are  Loschmidt’s 
number  and  the  standard  density  for  (neutral)  argon,  p  °the  preshock  density  is 
calculated  directly  from  the  measured  quantities  p^  and1  T^  via  the  equation  of 
state.  “ 


There  are  in  principle  two  sources  of  information  regarding  the 
values  of  the  two  remaining  unknown  quantities,  p^  and  2  (or  for  any  pair  of 

the  three  quantities  p^,  x^  and  Ng  2  which  can  be  employed’to  determine  (p^  -l)Q). 


Firstly,  these  quantities  can  be  calculated  directly  from  the  shock 
Hugoniot  equations  for  conditions  of  Saba  equilibrium  corresponding  to  the  experi¬ 
mental  values  of  M  ,  p  and  T  using  electronic  excitation  partition  functions 
si  1 

which  include  an.  array  of  excited  state  terms  as  described  in  Section  4.3.  The 
accuracy  of  this  calculation  is  limited  principally  by  the  random  errors  in  the 
measurements  of  Mg  and  p^  which  have  been  estimated  to  be  +  1%  and  +  2-1/2 $ 

respectively  (see  Section  3*l)*  For  a  Mg  =  20  shock  into  argon  at  p^  =  1  mmHg 
and  Tn  =  298°K,  a  +  1%  error  in  M  generates  a  +  1$  error  in  p0,  a  +  2-\/2%  error 
in  x  and  a  +  b%  error  in  N  D.  Similarly  a  +  2-1/2%  error  in  p.  generates  a 


ko 


I 

/ 


At- 


It  would 


+1/8$  error  in  p  ,  a  +  1/8$  error  in  x?  and  a  +  2-1/2$  error  in  N 
therefore  seem  desirable  to  avoid  use  of  the  calculated  values  of  N 


therefore  seem  desirable  to  avoid  use  of  the  calculated  values  of  N  _  and  x- 
.  \  .  ,  .  e,2  2 

to  determine  (p^.  -  1)q  from  Eq.  4.4.1  due  to  the  relatively  large  errors  produced 

by  the  random  errors  incurred  in  measuring  and  M  .  In  addition  the  calculated 
value  N  ^  was  shown  in  Section  4.3  to  be  on  average  about  2$  below  the  experi¬ 
mental  vdlue  Ne®g  which  was  suggested  to  be  due  to  use  of  an  unrealistic  model  for 
terminating  the ’partition  function  summation. 

The  second  source  of  information  regarding  the  values  of  p  ,  x^  and 
N  0  is  provided  by  the  method  of  two  wavelength  interferometry  which  has  already 

been  employed  to  determine  N  0  (see  Sections  2.,3  and  4.3).  By  simultaneous 

solution  of  Eq.  2.2.3  for  A  and  A  ,  the  laser  fundamental  and  second  harmonic 
^  y  x7 

wavelengths  respectively,  expressions  for  the  experimental  values  of  p  ,  x^  and 
0  are  obtained  directly  in  terms  qf  the  refractivities  for  Wie  component 
species  at  STP  and  the  experimentally  measured  fringe  shifts  S^  and  S^  .  Un¬ 
fortunately  the  expressions  for  p0  and  x0  are  dependent  on  the  'as  yet'  unknown 
X  f.  d 

values  of  (p^.  -  l)^  whereas  that  for  Ng  2  (Eq.  2.3.2)  is  to  a  first  approximation 
independent  of  the  ion  refractivity  as  discussed  in  Section  2.3. 

In  view  of  these  foregoing  comments  and  those  already  made  in  Section 
4.3-regarding  N  _,  it  was  decided  to  use  the  experimental  values  of  N  D  in  the 

expression  for  (  (p.^.  -l)Q  together  with  the  calculated  values  of  as  this 

quantity  has  been  shown  to  be  relatively  insensitive  to  errors  in  M  and  p1 . 

The  mean  values  for  the  refractivity  and  polarizability  of  singly 

ionized  argon,  calculated  from  the  experimental  data  via  Eq.  4.4.1  are  given  in 

Table  3  for^ wavelengths  of  3471$  and  5300$.  Also  shown  are  the  corresponding 

values  of  KT,  the  specific  refractivity  or  Gladstone-Dale  constant.  K”,  the 
X  s 

specific  refractivity  for  species  s  at  wave-length  A  is  related  tfo  the  refractivity 

(p£  -l)  through  the  relationship 


-  <*s  -  D 


(4.4.2) 


As  the  mass  densities  for  neutral  and  singly  ionized  argon  at  STP  are  equal  if 

one  neglects  the  difference  corresponding  to  the  mass  of  one  electron,  then 

throygh  use  of  the  equalities  represented  by  Eq.  1.1.1  and  Eq.  4.4.2,  we  have 

that  ,  ,  , 

/A  -\  t,A  /^A\ 


K  -  a. 

<^A  - 


( ah 

'  A7o 


As  this  ratio  is  frequently  employed  in  the  interferometry  of  argon  plasmas,  the 
values  for  A  =  3471$  and  5300A  are  also  given  in  Table  3.  The, individual  values 
for  the  ratio  of  the  ion  refractivity  to  that  for' the  atom,  (pj-l)o/(p^-l)o,  are 

presented  in  histogram  form  in  Fig.  19  and  20  for  A  =  3471$  and  5300$  respectively 
and  for  cell  widths  of  0.1,  together  with  the  corresponding  Gaussian  probability 
distribution  curves  (Ref.  77).-  The  histogram  for  the  23  measurements  of 
(p!J  -1)  /(p.  -l)  for  A  =  3471a,  exhibits  a  reasonable  correspondence  with  the 

Gaussian  curve  in  view  of  the  ^Limited  number  of  samples.  The  histogram  for 
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A  =  5300$  is  however  somewhat  more  erratic  reflecting  the  fact  that  only  seven 
results  were  available. 

The  large  spread  of  the  data  in  both  of  these  histograms  is  due 
primarily  to  the  fact  that  the  fringe  shift  contribution  for  the  ion  is  generally 
a  small  fraction  of  the  total  due  to  all  of  the  contributing  species.  For 
A  =  3*171$  and  shock  conditions  of  M  =  20,  p  =  1  mmHg  and  T.,  =  298  K  and  taking 
(p^  -  l)o/(n^  -  l)^  =  O.65,  the  ionscontributes  about  8 $  to  the  sum  of  the  moduli 

of  the  fringe  shifts  for  all  of  the  contributing  species.  Clearly  for  longer 
wavelengths,  this  ion  fringe  shift  contribution  will  become  an  even  smaller  fraction 
of,  this  total However  the  averages  of  the  individual  experimental  errors  for 
(p  -  l)  /  (p .  -  l)  calculated  using  the  estimated  measurement  errors  for  M  ,  p. , 

T^  and  §s  given  in  Section  3*3>  were  found  to  be  about  2.6  and  *j.7  times  the 

standard  deviation  at  3471$  and  5300$  respectively,  suggesting  that  the  random 
error  limits  given  in  Section  3*1*4  are  overly  pessimistic. 

Experimental  values  of  (p^  -  l)Q/(p^  -  l)Q  were  also  determined  for 

the  laser  fundamental  wavelengths  A  =  6943$  and  10,600°  A  and  were  found  to  be 
almost  exactly  equal  to  values  for  the  corresponding  laser  second  harmonic  wave¬ 
lengths.  This  degree  of  correspondence  was  regarded  with  some  suspicion  particularly 
as  the  fringe  shifts  for  these  relatively  long  wavelengths  are  due  principally  to 
the  contribution  from  the  free  electrons.  Investigation  cf  this  behaviour  showed 
that  the  apparent  agreement  was  a  direct  result  of  using  the  experimental  values 
of  2  to  obtain  the  ion  refractivity  from  Eq.  4.4.1.  In  deriving  the  expression 

for  N5  0  in  Section  2.3  in  terms  of  the  experimental  fringe  shifts,  it  was  assumed 
that  the  ratio  (p^  -  l)o/(pA  -  l)Q  was  independent  of  wavelength.  Although  large 

(possible)  deviations  from  this  assumption  were  shown  to  have  a  negligible  effect 

on  the  experimental  values  of  N  _  obtained  via  Eq.  2.3*2,  it  was  found  that  the 

6  ^  ^ 

effect  of  this  assumption  was  being  carried  through  into  the  experimental  values  of 
(Pj  -  4)q/(pa  “  l)Q  due  to  the  high  sensitivity  of  the  measurements  to  small  changes 

in  electron  density.  A  test  calculation  for  shock  conditions  of  M  =  20.,  p^  =  1 
mmHg  and  T,  =  298  K  indicated  that  a  real  +  %  difference  between  the  values  of 
(p  -  l)  /(p  -  l)  for  A  =  3471$  and  6943$  would  result  in  a  +  l°jo  error  in  this 

same  ratio  at  3471$  but  a  +  8$  error  at  6943$  if  this  difference  was  left  un¬ 
accounted  for  in  the  determination  of  the  experimental  values  of  If  „  through 

e,a 

Eq.  2.3.2  instead  of  the  more  exact  formulation  given  by  Eq.  2.3*1*  As  these 
percentages  indicate,  this  assumption  tends  to  bring  the  calculated  values  of 

(p^  -  l) J (p^  -■  i)^  for  the  two  wavelengths  Ax  and  Ay  into  close  agreement  while 

increasing  the  absolute  error  of  these  calculated  values  with  respect^to  the  real 
value.  It  is  therefore  obvious  that  the  assumption  that  (pj  -  1)q/-(pa  -  l)Q  is 

a  constant  in  the  derivation  of  2  from  the  experimental  data,  can  have  a 
significant  effect  on  the  ion  refractivity  at  the  (longer)  laser  fundamental 
wavelength  such  that  these  two  values  tend  to  asymptote  to  a  common  value  close 
to  the  real  value  for  the  wavelength  of  the  laser  second  harmonic.  The  ion  re-r 
fractivity  results  for  the  laser  fundamental  wavelengths  have  therefore  been  re¬ 
jected  due  ill  part  to  the  low  experimental  sensitivity  attainable  under  these 
conditions  and  in  part  to  the  tendency  of  these  results  to  take  on  values  close 
to  those  for  the  (shorter)  second  harmonic  wavelengths. 
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However  the  mean  experimental  values  for  the  ratio  of  the  ion  re¬ 
fractivity  to  that  for  the  neutral  atom  obtained  from  the  independent  interfero¬ 
metric  measurements  at  3471A  and  5300A,  are  in  close  agreement  (see  Table  3)  i.e., 
X  A 

(Pj  -  l)o/(uA  -  :)  =  0.65.  This  result  is  not  surprising  for,  as .discussed  in 

Section  2.1.,  the  dispersion  for  singly  ionized  argon  in  the  visible  and  near 
ultra-violet  parts  of  the  spectrum  should  be  similar  to  that  for  neutral  argon  as 
the  resonance  transitions  lie  in  the  yacuum  ultra-violet  region.  This  result 
therefore  validates  the  assumption  used  to  obtain  Ke,g  from  Eq.  2.3.2  that  the 

ratio  (pj  -  l)o/(pA  -  l)Q  is  essentially  independent* of  wavelength,  which  in  turn 


justifies  the  use  of  the  experimental  values  of 
from  Eq.  4.4.1. 


e,2 


to  derive  the  ion  refr activity 


It  is  now  relevant  to  discuss  the  possible  effects  of  the  excited 
electronic  states  on  the  argon  ion  fefractivity  results.  In  Section  2.1  it  was 
shown  that  the  excited  states  for  neutral  argon  make  contributions  of  approximately 
-1/2$  and  -1$  to  the  total  neutral  atom  refractivity  at  3471a  aqd  5300/frespecti- 
vely  for  an  argon  plasma  in  LTE  at  13 3?00°K  and  1  atmosphere  pressure  (see  Fig'.  I 
and  Table  2.).  As  these  contributions  were  not  included  in  the*  data  analysis  for 

(p.j.  -  l)Q,  they  will  represent  corresponding  corrections  to  (p^  -  1)q  of  the  order 

of  +3/2$  and  +  1$  respectively.  However  as  these  corrections  were  found  to  be- 
relatively  small  fractions  of  the  total  error  limits  given  for  (p^  -  l)  ,  they 
were  not  incorporated  into  the  final  refractivity  values  given  in  Table  3* 

In  addition,  as  already  discussed  in  Section  4.3?  the  consistency 
of  the  experimental  Ng  „  results  obtained  using  the  dual  frequency  ruby  and  neo¬ 
dymium  glass  lasers  indicated  that  the  effect  of  the  excited  states  was  not 
apparent  at  the  laser  fundamental  wavelengths.  The  only  significant  calculated 
excited  state  contribution  was  found  to  occur  at  the  ruby  laser  wavelength, 

A  =  6943X.  For  an  argon  plasma  in  LTE  (Ref.  75)  at  13,50O°K  and  1  atmosphere 
pressure,  the  excited  states  for  neutral  argon  were  calculated  to  produce  a  6.6$ 
reduction  in  the  total  species  refractivity  at  6943A.  As  shown  in  Section  4,3 
this  would  reduce  the  experimentally  measured  electron  densities  by  about  1$. 

Hence  as  the  same  calculation  indicated  a  negligible  contribution  at  the  neodymium 
glass  laser  wavelength  A  =  10,600X  and  us  the  average  values  for  N®  g  obtained 

using  the  two  lasers  differ  by  about  0.3$  in  a  direction  inconsistent  with  the 
calculated  excited  state  contribution,  it  would  seem  that  the  excited  states  have 
not  influenced  the  experimental  electron  density  results  used  to  calculate  (pj  -  l)Q 

from  Eq.  4.4.1,  for  the  prevailing  plasma  conditions.  It  is  therefore  concluded 
that  for  an  argon  plasma  in  LTE  for  conditions  of  13,5Cp0K  and  pressures  of  the 
order  of  1  atmosphere,  the  excited  states  for  neutral  and  singly  ionized  argon 
make  a  negligible  contribution  to  the  overall  plasma  refractive  index  at  the 
laser  wavelengths  employed  in  these  experiments.  The  refractive  index  data  for 
singly  ionized  argon  given  in  Table  3  therefore  in  effect  represents  the  ground 
state  term  for  this  species. 

^  Finally  it  is  interesting  to  compare  the  experimentally  measured 
ratio  (pT  -  l)  /(p/  -  l)  =  O.65  with  the  only  available  theoretical  value  of 

0.72  calculated  for  infinite  wavelength  using  Slater  screening  constant  theory 
(Ref.  10).  However  this  correspondence  does  not  apf>Ty  to  the  individual  refrac- 
tivilies  for  as  indicated  in  Section  1,  this  theory  is  at  best  only  approximate 
particularly  as  the  absolute  calculated  values  for  both  neutral  and  singly  ionized 
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argon  are  about  2-1/2  times  larger  than  the.  corresponding  experimental  values, 
k.5  Conclusions 

Interferometric  measurements  made  in  the  visible  (5300&)  and  near 
ultraviolet  ( 3U71S)  regions  indicate  that  the  refractivity  and  electric  .dipole 
polarizability  for  singly  ionized  argon  has  a  value  of  about  6,65  of  the 
corresponding  value  for  neutral  argon, 

This  experimental  determination  shows  reasonable  correspondence 
with  the  calculated  ratio  for  the  static  polarizabilities  but  as  yet  there 
exists  no  rigorous  calculation  for  the  polarizability  of  singly  ionized  argon 
with  which  to  compare  this  experimental  value, 

The  constancy  of  the  measured  electron  densities  and  ion  refrac- 
tivities  obtained  at  different  wavelengths  suggest  that  the  effects  of  the 
excited  state  refractivities  have  not  been  observed  for  the  conditions  of  the 
present  experiments. 

Experimental  confirmation  for  the  validit^_of  the  excited  state 
refractivity  calculation  for  neutral  argon  could  be  obtained  by  performing  shock 
tube  argon  plasma  interferometry  close  to  the  broad  resonance  peaks  calculated 
to  occur  at  about  8200&  and  9100X. 
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Table  2.  Calculated  excited  state  refractivities  for  neutral  Argon 
as  a  funtion  of  temperature  T  and  wavelength  X 
fora  Plasma  pressure  p=  I  atm.  electronic  port  it  ion 
function  is  shown  in  parenthesis  .bejow  corres¬ 
ponding  temperature. 

Refractivity  for  neutral  (grouRd  state)  Argon  at 
STP ,  Do,  is  shown  on  last  row  for 
comparison . 


TABLE  3.  REFRACTIVITY,  POLARIZABILITY  8  SPECIFIC  REFR ACTIVITY 
RESULTS  FOR  SINGLY  IONISED  ARGON 


(I-  \r/) 


FIG.  X.  CALCULATED  NEUTRAL  ARGON  REFRACTIVITX  INCLUDING  EXCITED 
STATE  CONTRIBUTIONS  AS  A  FUNCTION-OP' WAVELENGTH  FOR 
PLASMA  CONDITIONS  OF  13,500°K  AND  X  ATMOSPHERE  PRESSURE 


y  -  o 

o 


FIG.  2.  SCHEMATIC  LAYOUT  OF  ELECTRONICS  FOR  SHOCK  VELOCITY  MEASUREMENT  AND  LASER 
OPERATION . 


(a)  RUN  #51,  Ms =20; 50 ,  p  =  1.228  mm  Hg,  ,  T,  =297. 0°K 
UPPER  TRACE  10 /x  sec /'cm,  I  V/cm  ' 

LOWER  TRACE  I  /x  sec/cm ,  I  v  /cm 

GAUGE  COVERED  WITH  TRANSLUCENT  LAYER  OF 

EPOXY  RESIN.  ’  *■ 


{ b)  RUN  #57,  Ms  =  22.40,  pf  =0.614  mm  Hg ,  T,  “ 295.3°  K 
UPPER  TRACE  10/xsec/cm,  1/2  v  /cm 
LOWER  TRACE  lyxsec/cm,  1/2  v /cm 
GAUGE  COVERED  WITH  OPAQUE  LAYER  OF  BLACK 
SILICONE  { RTV)  RUBBER. 

FIG.  3  EFFECT  OF  RADIATION  PRECURSOR  ON  PIEZO¬ 
ELECTRIC  SMOCK  DETECTOR 


FI6URE4.  GEOMETRY  OF  9  INCH  MACH-ZEHNDER  INTERFEROMETER 


FIG.  5(a).  INITIAL  ALIGNMENT  OF  COMPLETE  OPTICAL  SYSTEM  WITH  HELIUM-NEON  LASER. 


HIGH  PRESSURE 
MERCURY  LAMP 


FIG.  5(t): INITIAL  FRIN GE ' 'ADJUSTMENT " OF  INTERFEROMETER 
USING  MERCURY  VAPOUR  LAMP,. 


FIG.  5 ( c )>.  ILLUMINATION  OF  INTERFEROMETER  BY  MEANS 
OF  HELIUM-NEON  LASER’. 


FROSTED  GLASS  COLLIMATING 

TUNGSTEN  LAMP  OPTICS 


FIG.  5(d)-  DIFFUSE  ILLUMINATION  OF  INTERFEROMETER 
BY  TUNGSTEN  LAMP  FOR  OBTAINING  FOCUSED 
FRINGES . 


FIG;.  5  ( e  ) .  ALIGNMENT.  OF.  RUBY /NEODYMIUM  LASER  WITH 
'  -  INTERFEROMETER  SYSTEM  BY  MEANS  OF  HELIUM- 

NEON  LASER. 


MACH-ZEHNDER 

INTERFEROMETER  CAMERA 


FIG.  5{f)  SCHEMATIC  OF  DUAL  FREQUENCY  INTERFEROMETER  SYSTEM 


) 


SIMULTANEOUS  MONOCHROMATIC  INTER- 
FEROGRAMS  FOR  A  M$=  17.46  SHOCK  INTO 
ARGON  AT  Pj  =2.32  mm  Hg  WITHOUT  THE 
ADDITION  OF  HYDROGEN. 

(I)  X=347IA  ,  (ii)X=6943&. 

SHOCK  TRAVELLING  JO  LEFT.  VERTICAL 
APERTURE  7  INCHES,  AN  INCREASE  IN 
REFRACTIVE  INDEX  MOVES  THE  FRINGES 
UPWARD. 


FIS.  8  SIMULTANEOUS  MONOCHROMATIC  INTER - 

FEROGRAMS  FOR  A  Ms=  16.60  SHOCK  INTO 
NITROGEN  AT  p=2.26  mm  Hg.  (i)  X=347I°A 

(ii)  X=6943°A.  SHOCK  TRAVELLING  TO 
LEFT.  VERTICAL  APERTURE  IS  7 INCHES. 

AN  INCREASE  IN  REFRACTIVE  INDEX  MOVES 
THE  FRINGES  UPWARD. 


( 


FIG  8(a)  INTERFEROGRAM  FOR  Mg  =  11.6  SHOCK  INTO  OXYGEN  AT  p  =' 1*0.0: 

mmHg,  A  =  691*32.  SHOCK  FRONT  AT  EXTREME  LEFT  AND  MOVING 
TO  LEFT.  INCREASE  IN  REFRACTIVITY  MOVES  FRINGES  DOWNWARDS 


HEAVY  PARTICLE  T4PERATURE 


SBESaEBUtMWS 


•» 


FIG.  9.  HEAVY  PARTICLE  AND:  ELECTRON  TEMPERATURES  BEHIND  MACH  .NUMBER  1*8.0  INCIDENT 
SHOCK  FRONT  .IN  ARGON-,  p..  =•  -3=  mm  -Hg-,  '  T,-.  ■=•  296. 5°K  -(-FROM  REF;  60  )  . 


DISTANCE  BEHIND  SHOCK  FRONT  (rtsm) 


FIG.  10.  TEMPERATURE  AND 'DEGREE'-  'OF'TofflZ&TION  PRO¬ 
FILES  BEHIND  THE  SHOCK  FRONT  IN  HYDROGEN., 
p  =  3  mmHg.,  M  =  25  (FROM  REF;  U8  )*, 

X  8 


INTERFEROGRAMS  OF  SHOCKS  INTO  ARGON 
WITH  SMALL  PERCENTAGES  OF  HYDROGEN 
X =6943°  A. 

(i)  MS=!7.0I ,  p,  =2.S4  mmHg  +  0.033%  H2 

(ii)  MS=I7.  II ,  pj  “2.85  mm  Hg  +  0,2%  H2 
SHOCK  TRAVELLING  TO  LEFT.  VERTICAL 
APERTURE  IS  7  INCHES.  AN  INCREASE 
IN  REFRACTIVE  INDEX  MOVES  FRINGES 
UPWARD. 


0)  (ii) 


FIG.  12  SIMULTANEOUS  MONOCHROMATIC  INTERFEROr 
GRAMS  FOR  A  Ms= 17.15  SHOCK  INTO  ARGON 

+  0.4%  H2  AT  p,s2.88  mm  Hg  (i)  X=347IA 
(ii)  X=6943 K.  SHOCK  TRAVELLING  TO  LEFT. 
REFERENCE  WIRE  SPACING  IS  2  INCHES. 

AN  INCREASE  iN  REFRACTIVE  INDEX  MOVES 
THE  FRINGES  UPWARD. 


EXPERIMENTAL  VALUE  FOR  PURE  ARGON  (REF-  58) 
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SHOCK  MACH  NUMBER  M, 


FIG,  3.1*1  EFFECT  OF  ADDITION  OF  O.tyf  HYDROGEN  ON  FROZEN  SHOCK 
PROPERTIES  FOR  SHOCK  MOVING  INTO  ARGON  SHOWN • AS 
PERCENTAGE  DEVIATIONS  FROM  THE  VALUES 'FOR  PURE  ARGON 
HYDROGEN  ASSUMED.  TO  /BE  'COMPLETELY  DISSOCIATED. 


FIG.  15  SIMULTANEOUS  MONOCHROMATIC  INTER FEROGRAMS 
FOR  A  Mg=  18.14  SHOCK  INTO  ARGON  +  0.4%  H2 
AT  p(= 2.45  mm  Hg  (i)X=5300°A  (ii)X  =I0600°A 
SHOCK  TRAVELLING  TO  LEFT.  REFERENCE  WIRE 
SPACING.  IS  -2  INCHES.  .  AN.  INCREASE  IN  REFRAC¬ 
TIVE  INDEX  MOVES  THE  FRINGES  UPWARD. 


DEGREE  OF  IONISATION 


COMPARISON  BETWEEN.  EXPERIMENTAL  ELECTRON  DENSITY  Ne  AND  Ns  CALCULATED  FOR 
SAHA  EQUILIBRIUM  BEHIND  INCIDENT  SHOCK  INTO  ARGON  USING  FULL  ELECTRON  EX¬ 
CITATION  PARTITION  FUNCTION..  ALSO  COMPARISON.  BETWEEN  NeF  and  NeG  CALCULATED 
AS  FOR  NeF  BUT  USING  ONLY  GROUND  STATE'  TERM  OF  PARTITION  FUNCTION. 


T,  «  298  *K 


SHOCK  .MACH  NUMBER  .  .M,  ....A 


FIG.  18.  ELECTRON  DENSITY  N  F  CALCULATED  FOR; 'SAHA 

EQUILIBRIUM  BEHIND  INCIDENT  SHOCK  INTO  ARGON 
USING  FULL  ELECTRONIC  EXCITATION  PARTITION 
FUNCTION. 


NUMBER  OF  MEASUREMENTS 


FIG.  19:  HISTOGRAM  OF  -  EXPERIMENTAL  (Hj-l)0/(^-l)0 

RESCUES  FOR  X®  3411  S'.'  MEAH  VAUJE  ®  0.650. 
TOUAIiJ  HUMBER  OF  'MEASUREMEHTS'  H  =  23; 
•STAHMED  DEVIATION  •0,l60.  SIAHDARD  • 
ERROR  OF  MEAH  (cr/\®)  *  0.033;' 


NUMBER  OF  MEASUREMENTS 


FIG.20:  HISTOGRAM  OF  EXPERIMENT Ai  0/,(^"X)0  BESETS  TOR 

Tv  5300  ;S.  MEAN  VALUE  '=0;647..  -TOTAL.  NUMBER  OF 
MEASUREMENTS  N  =  7.  STANDARD.  DEVIATION.  tr  =  0.107. 
STANDARD  ERROR  0?  .MEAN  (a/Ja)  =,;0.(Ali 


APPENDIX  A:  PROCEDURE  FOR  THE  ANALYSIS  OF  THE  INTERFEROGRAMS 

This  appendix  describes  the  procedure  -for  determining  the  equilibrium 
fringe  shift  between  pre-shock  state  1  and  post  shock  (equilibrium)  state  ,2  with 
the  shockwave  moving  into  the  region  of  state  1  in  a  direction  parallel  to  the 
fringe  field  so  that  the  shock  front  lies  perpendicular  to  this  same  fringe  field 
(see  Fig.  A-l(i.',)). 

For  an  ideal  interferometer,  devoid  of  inhomogeneities  due  to  optical 
defects,  vibrations  and  air  currents,  the  measurement  of  the  fringe  shift  is 
straightforward.  This  is  achieved  by  extrapolating  any  given  monochromatic  fringe 
in  region  1  into  the  desired  equilibrium  state  in  region  2  and  measuring  the  total 
shift  in  terms  of  the  fringe  spacing. 

In  the  present  system,  due  to.  imperfections  in  the  optical  compo¬ 
nents  such  that  the  fringes  connecting  regions  1  and  2  are  not  perfectly  straight 
and  horizontal,  it  is  necessary  to  take  an  interferogram  of  the  optical  field  in 
the  absence  of  the  shock  induced  fringe  shift.  This  therefore  provides  a  map  of 
the  optical  inhomogeneities  present  in  the  interferometer  which  can  then  be  sub¬ 
tracted  from  the  fringe  shift  measurements  obtained  from  the  shock  flow  inter-, 
ferograms  so  as  to  give  the  fringe  shift  due  to  the  shock  flow  alone.  It  is 
therefore  important  that  neither  the  flow  nor  the  no-flow  interferograms  contain 
any  fringe  movements  or  distortions  due  to  mechanical  vibrations  or  air  currents 
as  correction  for  these  shifts  is  beyond  the  control  of  the  present  experiment. 

The  interferometer  is  therefore-  sealed  in  an  -airtight  enclosure  and  care  is.  taken 
to  ensure  that  all  sources  of  low  frequency  vibrations  are  inoperative  during 
the  period  of  any  one  experiment.  ...  '  . 

However,  due  to  the  finite  time  which  exists  between  the  taking 
of  the  no-flow  and  flew  interferograms  minutes),  it  is,  possible  for  the  1 
interferometer  setting  to  change  due  to  thermal  or  mechanical  disturbances.  Thi3 
can  result  in  small  changes  in  fringe  spacing,,  position  pr .orientation;  thereby 
constituting  a  source  of  error  which  must  be  accounted  for  iaaddition  that  that 
due  to  the  inhomogeneities  of  the-  optical  system.  .  .  , 

Figure  A-l  shows  in,  schematic  form  the  flow  ,(i)  and  the  no-flow 
(ii)  interferograms  typical  of  those  encountered  in  this  work  where  three  dark 
monochromatic  fringes  are  represented  by  continuous  lines.  The  three  fringes 
in  the  no-flow  interferogram.  which  are  not  necessarily  the  same  three  in  the 
flow  interferogram  are  intentionally  given  an  exaggeratfedunon- ideal  appearance 
for  the  sake  of  this  discussion. 

The  broken  lines  represent  the  horizontal  and  vertical  reference 
wires  referred  to  in  Section  3.2  where  the  vertical  wire  is  located  so  as  to 
overlie  the  pre-shock  zone  (state  l)  at  a  point  close  to  the  shock  front. 

Point  X'  on  the  flow  interferogram  is  chosen  go  lie  at  any  desired  equilibrium 
(or  non- equilibrium)  state  2.  Point  Y',  on  the  no.-flow  interferogram  is  then 
chosen  such  that  YY'  =  XX *  on  the  assumption  that  the  two  interferograms  are 
reproduced  on  the  same  scale  otherwise  a  scaling -factor  is  used.., 

The  exercise  therefore  involves  the  measurement  of  the  fringe 
fractions  a,  b,  c?  d,  k,  1,  m  and  n  formed  on -either  side  of  the  reference  points 
X,  X’,  Y  and  Y' .  The  purpose  of  measuring  these  fringe  fractions,  rather  than 
following  a  given  fringe  through  the  sfeock  front  is  that  the  measurements  are 
localized  in  a  given  region  of  the  aperture  with  respect  to  the  given  fixed 


A1 


reference  points  X,  X' ,  Y  end  Y' .  Therefore,  provided  the  fringe  spacing  is  small 
in  relation  to  any  region  of  fringe  distortion,  it  should  be  possible  to  correct 
for  these  inhomogeneities  in  spite  of  small  differences  of  fringe  setting  which 
occur  between  the,  fringes  in  i egion  1  of  the  flow  interferogram  to  those  in  region 
1  of  the  no-flow  interferogram. 

The  shock  flow  fringe  shift  S^2(F)  between  pre- shock  state  1  and  shock 
state  2  can  be  given  as 


or  alternatively  as 


S12(F)  -  A  +  (c+d)  -  (a+bj 


s£2(f)  =  A  + 


d 

(a+c) 


(A-l) 


where  the  integral  fringe  shift  A  is  the  minimum  number  of  dark  fringes  one  crosses 
in  passing  from  X  to  X' ,  It  should  be  noted  that  both  equations  in  Eq.  A-l  are 
exactly  equivalent  as  they  can  both  be  reduced  to  the  common  form 


3^(F)  =  A-l  + 


b.  c 

(a+b)  (c.+d) 


Going  from  X  to  X'  integral  fringe  shifts  upwards  are  arbitrarily  denoted  as  being 
positive  and  those  downwards  as  negative  but  where  the  fringe  fractions  ajb*c,d> 
etc  are  scalar  quantities.  For  the  trivial  example  shown  in  Fig.  A- I,  A  =  +  1, 

A  =  +  1  and  B  =  -1.  Clearly  if  a  given  fringe  should  recross  either  of  these 
icf'erence  lines  between  points  X  and  X1  or  Y  and  Y’  ,  then  it  should  not  be  included 
in  the  values  for  A  or  B. 

Ideally  no  change  In  fringe  setting  occurs  between  the. -taking  of  the 
no-flow  and  flow  inter ferograms1.  „  In  reality  this-  is  not  the  case  due  to  .the 
effects  of  vibrations,  air  currents  and  temperature  changes  such  that  a  j  k  and 
b  i.  With  a  vj-ew  to  improving  the  accuracy  of  the  measurements  of  3j2(F),  the 
two  values  for  S12(F),  given  by  Eq.  A-l  are  averaged  to  give 


sJ2(f)  =  A  +  \  j. 


:(A-2) 


In  a  similar  manner  it  is  now  possible  to1  obtain  the  corresponding 
vertical  fringe  shift  S^2(NF)  occurring  between  states  1  and  2  on  the  norflow 
interferogram  due  to  optical  defects  in  the  system.  The  equivalent  average  value 
for  S^9(NF)  is  given  by 


s£2(nf)  =  b  +  \- 


(A-3)’ 


The  true  total  fringe  shift  S^2  due  to  the  gas-dynamic  phenomenon  under,  investir 
gation  is  therefore  given  by 

4  =  4(p)  -  "#<-> 

which  by  substitution  from  Eq.  A-2  and  A-3  becomes 


A2 


S^2  =  A-B  +  1/2 


:>-4> 


where  the  denbiniuators  for  the  terms  within  the  brackets  represent  tiie  fringe  spacings 
about  reference  X,  X',  Y  and  Y'.'  Due  to  optical  defects  ih!  the  system,  these  fringe 
spacings  do  not  have  a  common  value  although  normal  circumstances  the  differences 

are  small.  *  "  ’•  '  ’  \  j ''  •*  '  1  - 

It  is  now  pertinent  to  discuss'  the  effect  of  large  optical  ihhomo- 
geneities  and  changes  of  fringe  setting  which  occdr  between  the  taking  of  the  no¬ 
flow  and  flow  inteiferograms  on  the  validity  of  this  analytical; 'procedure. -  £f 
the  fringe  spacirig  is  large  in  relation  to  any  given  region  of  optical  imperfection, 
then  it  is  possible  for  this  region  to  lie  between  two  dark  fringes  at  point 'X ’  on 
the  flow  interferogram  (say)  such  that  the  fringes  appear  to  be  straight  and  parallel. 
However,  this  same  region  may  lie  at  some  intermediate  fringe  position  at  point  Y1 
on  the  no-flow  interferogram  such  that  the,  local  fringe  profile  now  exhibits  a  non¬ 
ideal  shape.  It  is  therefore  clear  that  the  no-flow  interferogram  will  introduce 
a  fringe  shift  correction  for  an  optical  defect  which  is  not  included  in  the  data 
from  the  flow  interferogram.  The  effect  of  this  anomaly  can  be  minimized  by  taking 
one  or  both  of  two  precautions.  Firstly  by  ensuring  that  the  fringe  spacing  is  . 
much  smaller  than  any  region  of  optical  inhomogeneity,  the  fringe  field  effectively 
becomes  a  continuous  map  of  the  transmitted  optical  wavefront  with  all  defects  being 
clearly  defined.  This  expedient  is  of  course  subject  to  the  proviso,  that  the  fringe 
spacing  be  much  larger  than  either  the  photographic  emulsion  or  diffuser  generated 
grain  size  or  the  reference  wire  diameter. 

The  other  means  by  which  this  difference  between  the  two  interferograms 
can  be  minimized  is  by  taking  an  average .value  for  the  fringe  spacing  about  points 
X  and  Y  and  similarly  for  points  X1  and  Y’ .  These  average  spacings  are  (c  +  d  + 
m  +  n)/2  and  (b  +  a  +  k  +  i)/2,  which  when  substituted  for  the  original,  values  used 
in  Eq.  A-2  and  A-3,  provide  a  new  expression  for  the  total  fringe  shift  given  by 


4  -  A-B  -  )  <A-5> 

However  except  for  relatively  large  differences  between  (c+d)  and  (m+n)  and  between 
(b+a)  and(l+k),  the  value  of  S^2  given  by  Eq.  A- 5  is  not  likely  to  differ  from  the 

value  given  by  Eq.  A-4  by  more  than  a  fraction  of  one  percent. 

The  effect  of  a  change  of  fringe  setting  between  the  no-flow  and  flow 
interferograms  on  the  value  of  can  be  considered  in  three  parts.  Firstly  a 

change  in  vertical  fringe  position  will  only  affect-  the  fringes  for  preshock  state 
1  as  a  change  in  fringe  position  between  the  no-flow  and  flow  interferograms  for 
post-shock  state  2  already  exists.  Clearly,  the  effect  of  this  movement  is  subject 
to  the  comments  already  made  regarding  this  point  for  region  2 .  Secondly,  a  change 
in  fringe  spacing  (  <1/20  fringe  say)  will  not  have  a  noticeable  effect  on  S^V,  as 
calculated  from  either  Eq.  A-4  or  A-5  provided  the  fringe  spacing  is  sufficiently 
fine  to  resolve  all  the  optical  imperfect  ions  in  the  system.  This  is  due-  to  the 
fact  that  the  ratio  of  the  fringe  fractions  to  the  fringe  spacings  remains  constant. 

Finally,  the  effect  of  a  change  in  fringe  orientation  away  from  the 
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horizontal  setting  occurring  between  the  no-flow  and  flow  interferograrus  cannot  be 
dismissed  as  easily  as  for  the  effect,  of  the  other  changes  in  fringe  adjustment 
already  discussed.  Besides  generating  an  apparent  change  in  fringe  spacing  as 
viewed  in  the  vertical  direction,  the  measured  fringe  shift  for  the  flow  interfero- 
grams  will  clearly  depend  on  the  length  of  XX' .  Precautions  should  therefore  be 
taken  to  minimize  this  change  in  fringe  orientation  from  the  horizontal  direction 
by  isolating  the  interferometer  from  all  external  influences  and  by  reducing  the 
time  -period  between  the  taking  of  the  no-flow  and  flow  interferograrus .  In  addition 
the  effect  of  this  change  on  can  be  minimized  by  making  XX' as  small  as  possible 
or,  in  the  presence  of  large  changes  in  fringe,  orientation  of  the  order  of  5°  to 
10°,  by  making  the  fringe  shift  measurements  in.  a  direction  perpendicular'  to  the 
given  fringe  field  instead  of  in  a.  direction  perpendicular  to,  the  horizontal  ref¬ 
erence  wire,,  so  that  the  fringe  fractions  a,  b,  c  and.  d  maintain  ah  approximately 
correct  numerical  relationship  to  the  values  k,  1,  m  and  n.  for  the  no-flow  interf ero- 
gram. 


